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1.  ABSTRACT 


This  project  is  inserted  within  the  scope  of  the  technology  of  morphing  aircrafts.  The  research  team 
investigated,  in  a  multidisciplinary  framework,  the  design  methodology  of  variable  camber  morph¬ 
ing  wings  based  on  the  use  of  active  materials,  namely  piezoelectric  materials  and  shape  memory 
alloys.  The  research  work  has  been  carried  out  in  partnership  with  Prof.  Daniel  J.  Inman,  from  the 
Department  of  Aerospace  Engineering  of  the  University  of  Michigan.  After  preliminary  assessment 
of  various  innovative  candidate  morphing  strategies,  the  research  evolved  to  the  concept  of  active 
meta-structures  based  on  shape  memory  alloys,  which  is  currently  in  development  in  the  context  of 
a  Ph.D.  thesis  jointly  advised  by  Profs.  Domingos  A.  Rade  and  Daniel  J.  Inman. 

The  main  accomplishments  made  are  the  following: 

1 .  development  of  a  methodology  for  the  multidisciplinary  design  of  variable  camber  morphing 
wings  based  on  finite  element  simulations,  encompassing  the  modeling  of  the  thermo-mechanical 
behavior  of  shape  memory  alloys  and  aerodynamic  transient  analysis.  This  methodology  enables 
to  evaluate  the  evolution  of  the  aerodynamic  forces  and  moment  during  the  morphing  transition. 

2.  conceptual  design  and  tests  of  variable  camber  aerodynamic  profiles  combining  shape 
memory  alloy  wires  and  permanent  magnets  to  achieve  multi-stable  configuration  for  morphing 
purposes. 

In  the  sequence,  the  research  methodology  and  obtained  results  are  presented  and  discussed. 
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2.  OBJECTIVES 


The  main  objective  of  the  original  proposal  remained  unchanged,  consisting  in  the  investigation  of 
innovative  contributions  to  the  technology  of  morphing  aircrafts,  exploring  smart  materials,  aiming 
at  achieving  more  efficient  morphing  capability  in  terms  of  control  authority  and  energy  consump¬ 
tion. 

Other  specific  objectives  established  for  the  project  are: 

•  The  acquisition,  production  and  difusion  of  knowledge  on  a  theme  that  holds  great  scien¬ 
tific  and  technological  innovation  potentials; 

•  The  strengthening  of  the  scientific  cooperation  between  Brazilian  and  American  research 
teams; 

•  The  production  of  high-impact  publications. 

RESEARCH  TEAM 


The  research  team  is  composed  by  the  members  listed  below.  As  can  be  seen,  the  team  is  com¬ 
posed  by  researchers  of  various  academic  levels,  ranging  from  undergraduate  students  to  profes¬ 
sors. 


Name 


Position* 


Domingos  Alves  Rade 
Francisco  Jose  de  Souza 
Willian  Mota  Baldoino 
Thiago  de  Paula  Sales 
Leonardo  Sanches 
Vergilio  Torezan  Silingardi  Del  Claro 
Luigi  Albieri  Antonialli 
Bruno  Barcelo 

Alexandre  Felipe  Medina  Correa 


Professor,  PI 

Associate  Professor,  Co-PI 
M.Sc.,  Researcher 
M.Sc.,  Doctorate  Student 
Ph.D.,  Post-Doctoral  Researcher 
Undergraduate  Student 
Undergraduate  Student 
Undergraduate  Student 
Undergraduate  Student 


3.  PROJECT  ACTIVITIES  AND  RESULTS 


The  research  activities  have  been  conducted  with  some  modifications  with  respect  to  the  origi¬ 
nal  proposal,  according  to  the  following  steps: 
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Literature  review 


A  comprehensive  literature  review  was  carried  out  to  characterize  the  state-  of-the-art  related  to 
the  morphing  aircraft  technology  and  smart  structures.  Emphasis  has  been  put  on  the  following 
aspects:  i)  strategies  employed  for  adaptive  structural  systems,  their  advantages  and  disad¬ 
vantages;  ii)  multidisciplinary  numerical  modeling  techniques  and  experimental  methodologies 
needed  for  the  design  and  assessment  of  the  performance  of  morphing  wings.  The  literature  re¬ 
view  as  accomplished  by  the  study  of  a  large  number  of  conference  and  journal  papers,  and  tech¬ 
nical  reports  by  the  members. 

Study  mission  at  the  University  of  Michigan 

Four  members  of  the  research  team  of  this  project  have  made  a  study  mission  in  the  USA  :  Prof. 
Domingos  Alves  Rade,  Dr.  Francisco  Jose  de  Souza,  Mr.  Willian  Mota  Baldoino,  and  Mr.  Thiago 
de  Paula  Sales.  Professor  Valder  Steffen  Jr.,  a  faculty  of  the  School  of  Mechanical  Engineering  of 
the  Federal  University  of  Uberlandia,  has  joined  the  group  in  this  study  mission,  but  his  travel  has 
not  been  funded  by  AFOSR. 

The  study  mission  has  comprised  the  visit  to  three  research  groups  from  three  American  Universi¬ 
ties,  which  develop  research  devoted  to  morphing  aircraft,  namely: 

•  University  of  Michigan  at  Ann  Arbor,  Ml  (April  28  and  29,  2013).  The  team  was  received 
by  Prof.  Daniel  J.  Inman,  Dean  of  Department  of  the  Aerospace  Engineering,  as  well  as 
by  Prof.  Carlos  Cesnik.  Meetings  have  been  arranged  with  Amy  Conger  and  Alec  D.  Gal- 
limore,  who  were,  at  the  time  of  the  visit,  Director  of  international  programs  of  the  College 
of  Engineering,  and  Associate  Dean  for  research  and  graduate  education,  respectively. 
The  facilities  of  the  College  of  Engineering  of  the  University  of  Michigan  were  visited  by 
group.  Lectures  were  given  by  Prof.  Domingos  A.  Rade,  Prof.  Valder  Steffen  Jr.,  and 
Dr.Francisco  Jose  de  Souza  to  the  professors  that  welcomed  the  team,  so  that  common 
research  interests  could  be  better  established; 

•  Virginia  Polytechnic  Institute  and  State  University  at  Blacksburg,  VA  (May  1,  2013).  The 
team  was  received  by  Prof.  Pablo  A.  Tarazaga.  Formal  meetings  took  place  with  the  Head 
of  Department  of  Mechanical  Engineering,  Prof.  Robert  G.  Parker,  as  well  as  with  Prof. 
Glenda  R.  Scales,  the  Associate  Dean  for  International  Programs  and  Information  Tech¬ 
nology  by  the  time  of  the  visit.  Several  laboratories  have  been  visited,  including  the  Center 
for  Intelligent  Material  Systems  and  Structures  (CIMMS).  Lectures  were  given  by  Prof. 
Domingos  A.  Rade,  Prof.  Valder  Steffen  Jr.,  and  Dr.  Francisco  Jose  de  Souza  during  the 
visit; 
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•  Florida  State  University  and  the  Florida  A&M  University,  at  Tallahassee,  FL  (May  3,  2013). 
The  welcoming  team  of  the  FCAAP  comprised  Prof.  Rajan  Kumar,  Prof.  Kunihiko  Taira, 
Prof.  Chian  Shih,  Prof.  Bruce  R.  Locke,  and  Prof.  William  S.  Oates.  Professors  of  the  visit¬ 
ing  group  and  from  FCAAP  have  shared  their  common  interests  through  lectures,  as  well 
as  during  the  walk  through  some  of  the  research  laboratories  belonging  to  the  FCAAP. 

As  an  immediate  outcome  of  the  visit,  it  was  agreed  that  one  of  the  members  of  the  visiting  group, 
Thiago  de  Paula  Sales,  who  is  advised  by  Prof.  Domingos  A.  Rade,  will  perform  a  research  ex¬ 
change  program  at  the  University  of  Michigan  at  Ann  Arbor.  Starting  from  January  2015,  he  will  do 
part  of  his  PhD  studies  for  a  period  of  12  months  under  the  co-supervision  of  Prof.  Daniel  J.  Inman. 
At  the  time  this  report  is  written,  Mr.  Thiago  is  in  the  course  of  his  exchange  program,  funded  by 
the  Brazilian  Science  without  Borders  Program. 

Computer  simulation  of  variable  camber  morphing  wing 

During  this  stage,  the  team  got  acquainted  with  computational  software  for  the  multidisciplinary 
modeling  adaptive  airfoils.  Finite  element  software  ANSYS®,  which  was  made  available  by  the 
Federal  University  of  Uberlandia,  was  chosen  as  the  main  aerodynamic  and  structural  modeling. 
Specifically,  modules  of  the  ANSYS®  15.0  suite  Fluent®,  which  is  a  Computational  Fluid  Dynamics 
solver  that  incorporates  several  turbulence  models,  ANSYS  Mechanical®  and  ANSYS  Mechanical 
APDL®,  have  been  used  for  computer  simulations.  Moreover,  the  ANSYS  Workbench®  interface 
has  also  proven  useful. 

The  primary  goal  of  this  research  step  was  the  development  of  a  methodology  for  the  two- 
dimensional  transient  aerodynamic  modeling  of  variable  camber  airfoil,  which  enabled  to  evaluate 
the  evolution  of  aerodynamic  forces  and  moments  during  a  prescribed  shape  transition. 

The  above-mentioned  software  enables  the  creation  of  User  Defined  Functions  (UDFs),  which  are 
computer  codes  written  in  #C  language,  to  perform  various  user-defined  tasks.  One  of  the  features 
that  can  be  implemented  is  related  to  mesh  adaptation  as  required  through  the  code.  Thus,  a  rou¬ 
tine  based  on  UDFs  could  be  put  forward  to  simulate  transient  phenomena  related  to  adaptive, 
shape  changing  airfoils.  First,  a  mesh  must  be  generated  for  the  initial  geometry  for  the  problem. 
This  is  the  only  mesh  that  needs  to  be  input  to  the  software.  It  can  be  of  any  type  (hexagonal,  te¬ 
tragonal,  etc.)  and  generated  through  any  mesh  generation  software  as  long  as  it  format  is  ade¬ 
quately  recognized  by  the  solver.  Next,  the  written  UDFs  need  to  be  compiled  by  the  ANSYS  Flu¬ 
ent®  software,  a  procedure  that  may  entail  certain  difficulties  depending  on  the  operational  system 
and  the  lack  of  certain  functions  needed  by  the  compiler  in  the  default  root  of  the  system.  As  to  the 
UDF  #C  code,  it  can  be  developed  by  making  advantage  of  macros  listed  in  the  ANSYS  Fluent® 
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documentation  manual.  More  specifically,  the  UDFs  used  for  the  simulation  of  adaptive  geometry 
should  adopt  macros  that  define  a  dynamic  mesh,  i.e.,  a  mesh  that  can  change  throughout  time 
integration  of  field  variables.  The  most  general  approach  considers  grid  motion,  in  which  the  posi¬ 
tion  occupied  by  each  node  of  the  mesh  (or  of  a  region  of  it)  is  user-specified. 

Results  of  two  simulations  are  shown  here.  In  the  first,  a  NACA  0012  airfoil  has  its  shape  com¬ 
pletely  changed  to  that  of  a  Selig  1223  airfoil.  Both  profiles  are  shown  in  Fig.  1.  In  the  second,  a 
NACA  0012  airfoil  has  its  trailing  edge  rotated  about  the  80%  chord  position  by  20  degrees  down¬ 
wards.  Both  simulations  are  two-dimensional,  with  a  Reynolds  number  of  105  based  on  the  far-field 
speed,  and  on  a  chord  length  of  1  m.  The  initial  mesh  is  a  hexa  C-grid  that  was  created  using  AN- 
SYS  ICEM  CFD®;  it  is  fine  enough  to  capture  boundary  layer  effects  and  thus  better  predict  aero¬ 
dynamic  forces.  The  fluid  is  air,  and  the  k-co  SST  turbulence  model  is  adopted. 


a)  NACA  0012  b)  Selig  1223 


Figure  1  -  Illustration  of  two  airfoils  adopted  from  shape  transition  simulations 

The  positions  of  the  nodes  that  form  the  airfoil  surface  are  the  only  ones  that  are  updated  through 
UDFs.  The  remaining  nodes  of  the  mesh  are  automatically  moved  as  the  airfoil  surface  changes  its 
form  according  to  internal  algorithms  of  the  ANSYS  Fluent®  software.  These  can  be  chosen  by  the 
user,  and,  for  the  simulations  considered,  one  has  adopted  diffusion-based  smoothing,  a  technique 
that  can  maintain  the  mesh  geometry  unchanged  within  boundary  layer  regions  (thus  rendering 
better  quality  results). 

To  update  the  airfoil  surface  nodes,  a  time  interpolation  scheme  is  needed,  as  only  initial  and  final 
geometry  is  known  (specified).  This  is  not  a  UDF  restriction,  but  rather  an  example-related  issue.  If 
experimental  data  were  available,  these  could  be  directly  used  to  specify  node  positions  at  all 
times.  It  is  assumed  that: 


x(t)  =  xi+(xf-xi)j(t)\ 
y(t)  =  yi+(yf-yi)J(t), 

where  t  is  time,  (x,,y,),  ( xf,yf ),  and  (x,y)  denote  the  initial,  the  final,  and  the  time-dependent 
node  positions,  and  J(t)  is  the  interpolation  function  that  must  be  null  before  the  initial  time  for 
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shape  changing,  f, ,  and  unity  after  the  final  time  for  shape  changing,  tf .  Example  of  interpolation 
functions  are: 


•^linear  (  T)  T 

= T  -  is'ni-2'r‘Y- 

W  =  35r‘  -  84 r5  +  70r"  -20 r7, 

where  r  =  (f  -  f,.)/(ff  -  f() .  These  expressions  are  valid  for  0  <  r  <  1 ;  for  r  <  0 ,  J  =  0  ,  and  for  r  >  1 , 
J  =  1 .  Plots  of  the  different  interpolation  functions  are  depicted  in  Fig.  2. 


*5 


T 


Figure  2  -  Interpolation  functions  used  to  prescribe  the  motion  of  an  adaptive  airfoil. 


The  results  shown  in  Figs.  4  and  5  show  the  transient  nature  of  lift  and  drag  coefficients  (C;  and 
Cd ,  respectively)  of  both  considered  adaptive  airfoils.  Numerical  integration  was  made  with  an  im¬ 
plicit,  second-order  scheme.  The  responses  were  observed  for  a  relatively  long  time  period  to  as¬ 
sure  steady  state  equilibrium  after  geometric  modification. 


Time  [s]  Time  [s] 

Figure  3  -  Simulation  results  for  a  NACA0012  airfoil  evolving  to  a  Selig  1223  airfoil  considering 

different  interpolation  functions. 
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Figure  4  -  Simulation  results  for  a  NACA0012  that  has  its  trailing  edge  rotated  by  the  80%  chord 
location,  by  20  degrees  downwards,  using  polynomial  interpolation. 

As  one  can  see,  these  results  show  that  the  achievement  of  steady  state  response  depends  on 
the  adaptive  airfoil  shape  change.  When  geometrical  changes  are  relatively  mild,  as  when  the  trail¬ 
ing  edge  of  the  NACA  0012  is  simply  rotated  downwards,  the  lift  coefficient  curve  follows  the  inter¬ 
polation  function  behavior.  However,  when  strong  shape  changes  take  place,  as  in  the  case  in 
which  one  departs  from  a  NACA  0012  to  a  Selig  1223,  the  lift  coefficient  curve  displays  a  more 
complex  pattern,  leading  to  a  long  time  period  before  the  settlement  of  steady  state  conditions. 
This  observation  correlates  to  the  development  of  a  leading  edge  bubble  in  the  intrados  of  the 
Selig  1223  profile,  which  propagates  down  to  the  trailing  edge,  giving  rise  to  a  decrease  in  the  lift 
coefficient  value  after  geometric  shape  changes  have  ceased.  Hence,  one  can  see  that  direct  in¬ 
terpolation  of  aerodynamic  data  for  airfoils,  or  use  of  steady  state  data  to  model  transient  behavior 
might  fail.  Caution  should  be  taken  not  only  when  intense  shape  changes  are  observed,  but  also 
for  scenarios  in  which  pressure  gradients  assume  large  values.  (In  this  regard,  it  should  be  finally 
said  that  such  transient  effect  owned  to  pressure  gradients  might  be  shaped  to  give  superior  per¬ 
formance  of  airfoils,  through  constant  shape  modification  that  enables  controlled  separation  of 
bubbles  in  the  extrados  and  in  the  intrados  surfaces  -  thus  producing  a  net  upward  lift). 

Another  interesting  feature  that  comes  with  the  results  shown  is  that  too  drastic  shape  changes 
such  as  that  obtained  with  linear  interpolation  may  be  problematic.  If  the  responses  plots  are  truth¬ 
ful,  as  indicated  by  proper  convergence  of  residues,  impact  forces  would  be  directed  to  the  struc¬ 
ture  of  the  adaptive  airfoil.  This  phenomenon  could  lead  to  potential  problems  for  the  flight  control 
subsystem  of  an  actual  aircraft,  as  well  as  to  a  possible  weight  penalty  due  to  structural  reinforce¬ 
ments.  Thus,  by  the  presented  results,  actuation  for  morphing  airfoils  should  be  as  smooth  as  pos¬ 
sible.  Nevertheless,  further  assessment  is  still  to  be  made  regarding  this  conclusion  as  the  result¬ 
ing  discontinuities  in  the  C;  and  Cd  curves  might  be  owned  to  internal  algorithms  adopted  in  the 
ANSYS  Fluent®  software. 
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Figure  5  shows  graphically  a  typical  time  evolution  of  the  velocity  fields  around  the  shape  morphing 
airfoil. 


r  =  1,4  s  f=l.9s 

- - -  * 

_  ^ _ -  —  _  Jmm* 


f  =  1.5  s  _  t  =  2.0s 


Figure  6  -  Example  of  the  time  evolution  of  the  velocity  field  around  a  shape  morphing  airfoil. 


Construction  of  prototypes  and  wind  tunnel  tests. 

In  partnership  with  the  research  group  lead  by  Prof.  Carlos  J.  de  Araujo,  from  the  Federal  Universi¬ 
ty  of  Campina  Grande,  Brazil,  wind  tunnel  experiments  were  carried  out  on  prototypes  of  a  morph¬ 
ing  wing,  whose  articulated  ribs  could  have  their  camber  changed  by  the  activation  of  shape 
memory  springs  (Fig.  7). 


Figure  7  -  Articulated  wing  rib  actuated  by  shape  memory  springs. 
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The  tests  were  performed  at  the  Aeroelasticity  Laboratory  of  the  School  of  Mechanical  Engineering 
of  the  Federal  University  of  Uberlandia. 

The  main  goal  of  the  experiment  was  to  assess  the  performance  of  the  proposed  concept  in  terms 
of  the  time  necessary  for  transitioning  under  actual  aerodynamic  loads.  Images  shown  in  Fig.  8 
were  captured  during  the  wind  tunnel  tests,  which  have  validated  the  concept  but  showed  that,  as 
expected,  the  SMA  actuator  have  a  limited  to  counteract  the  aerodynamic  forces  that  arise  for 
higher  flow  speeds. 
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A  second  group  of  prototypes  has  been  developed  at  the  Federal  University  of  Uberlandia  as  part 
of  Mr.  Thiago  de  Paula  Sales’s  Ph.D.  thesis,  who  spend  one  year  at  the  Aerospace  Engineering 
Department  of  the  University  of  Michigan,  as  a  research  exchange  visitor,  under  the  supervision  of 
Dr.  Daniel  Inman. 

These  prototypes  have  been  developed  based  on  the  idea  of  combining  the  concept  of  meta¬ 
structures  with  shape  memory  alloys  and  permanent  magnets  to  achieve  multiple  stability  of 
morphing  airfoils.  Illustrations  of  these  prototypes  are  provided  in  Figs.  9  and  10. 


Fig.  9  -  Prototype  exploring  the  concepts  of  meta-structures  whose  unit-cells  are  driven  by  SMA 

wires  and  permanent  magnets. 


Fig.  9  -  Prototype  exploring  the  concepts  of  meta-structures  whose  unit-cells  are  driven  by  SMA 
wires  and  permanent  magnets  applied  to  a  morphing  airfoil. 
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FINAL  REMARKS 

The  research  team  evaluates  that  the  goals  originally  established  for  the  project  have  been  satis¬ 
factorily  attained.  Indeed,  the  phase  funded  by  AFOSR  has  enabled  the  Brazilian  research  group 
to  get  acquainted  with  advanced  aspects  of  the  morphing  wing  technology  and  develop  the  basic 
modeling  and  experimental  procedures.  More  importantly,  it  paved  the  way  for  a  far-reaching  re¬ 
search,  which  includes  the  development  of  a  Ph.D.  thesis  co-advised  by  profs.  D.A.  Rade  and  D.J. 
Inman.  From  this  thesis,  relevant  scientific  achievements  and  publications  are  expected  to  come. 
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Resumo:  Aeronaves  adaptativas  tem  recebido  grande  atengao  da  comunidade  internacional  pelos  beneffcios  que  podem 
propiciar.  Pode-se  citar  como  exemplo  a  possibilidade  de  redugao  no  consumo  de  combustivel,  devido  a  um  melhor  de- 
sempenho  aerodinamico  da  aeronave  em  seu  envelope  de  voo.  Apesar  de  varios  estudos  terem  considerado  o  projeto  de 
aeronaves  adaptativas,  levando  a  subsequente  fabricagao  e  realizagao  de  ensaios,  poucas  investigagoes  preocupam-se 
com  a  analise  do  periodo  transiente  no  qual  ocorre  reconfiguragao  geometrica  da  aeronave.  Neste  trabalho  sao  apre- 
sentadas  as  caracteristicas  aerodinamicas  (coeficientes  de  sustentagao,  arrasto  e  momento)  relacionadas  a  um  aerofolio 
adaptativo.  O  mesmo  assume  como  forma  inicial  o  perfil  NACA  0012,  e  evolui  para  o  perfil  Selig  1223.  Para  tanto,  simu- 
lagdes  bidimensionais  do  escoamento  ao  redor  do  aerofolio  sao  realizadas  com  auxilio  do  software  ANSYS  Fluent®  15.0. 
Uma  metodologia  de  parametrizagao  espacial  de  aerofolios  e  apresentada  e  utilizada  para  implementagao  das  mudan- 
gas  geometricas  envolvidas.  Os  resultados  obtidos  sao  comparados  com  outros  oriundos  de  simulagoes  que  consideram 
regime  permanente.  Os  mesmos  permitem  evidenciar  que  a  forma  com  que  se  da  a  alteragao  geometrica  do  aerofolio 
adaptativo  impacta  diretamente  sobre  os  coeficientes  aerodinamicos  computados,  podendo  acarretar  descontinuidades, 
e  eventuais  problemas  para  o  subsistema  de  controle  de  uma  aeronave. 

Palavras-chave :  Aeronaves  Adaptativas,  Fluidodinamica  Computacional. 


1.  INTRODUCE  AO 

Aeronaves  adaptativas,  ou  morphing  aircraft ,  em  ingles,  sao  caracterizadas  por  alteragoes  controladas  de  geometria. 
Tais  alteragoes  sao  de  grande  amplitude,  ou  localizadas  e  suaves,  e  portanto  se  diferenciam  dos  dispositivos  atualmente 
empregados  com  fins  de  hiper- sustentagao  e  de  controle  de  aeronaves.  O  interesse  em  seu  estudo,  que  teve  imcio  ha 
mais  de  tres  decadas  ( e.g .  Gilbert,  1981),  e  atribuido:  a  possiveis  ganhos  em  desempenho,  o  que  pode  acarretar  menores 
impactos  ambientais  e  redugao  de  custos  operacionais;  e  tambem  a  possibilidade  de  concepgao  de  aeronaves  capazes  de 
atender  a  varios  propositos,  sobretudo  com  requisitos  conflitantes  (Bolonkin  e  Gilyard,  1999;  Weisshaar,  2013).  Para 
tanto,  as  ditas  aeronaves  adaptativas  geralmente  consideram  alteragoes  no  enflechamento,  na  envergadura,  na  torgao  e  no 
diedro  de  suas  asas,  e/ou  manipulagoes  no  arqueamento,  na  espessura  e  na  corda  de  segoes  que  a  compoem. 

Uma  variedade  de  artigos  cientificos  e  teses  ja  se  encontra  dispomvel  na  literatura  sobre  o  tema  (Sofia  et  al. ,  2010; 
Barbarino  et  al.,  2011;  Valasek,  2012).  Grandes  projetos  de  pesquisa  tambem  ja  foram  concluidos  (Wlezien  et  al.,  1998; 
McGowan  et  al.,  2002).  Atualmente,  grandes  empresas  tambem  tem  investido  esforgos  para  implementagao  de  tecnolo- 
gias  adaptativas,  no  sentido  descrito  anteriormente,  visando  aumento  da  eficiencia  energetica  de  aeronaves  (Butterworth- 
Hayes,  2014). 

Como  resultado  do  esforgo,  algumas  aeronaves  adaptativas  ja  foram  projetadas,  construidas  e  testadas  em  tuneis  de 
vento  e  em  voo  -  ainda  que  se  tratem,  em  sua  maioria,  de  provas  de  conceito,  e  sejam  de  pequena  escala  (Barbarino 
et  al.,  2011;  Gomez  e  Garcia,  2011).  Apesar  do  sucesso  de  trabalhos  ja  feitos,  o  comportamento  transiente  de  aeronaves 
adaptativas  e  relativamente  pouco  estudado  e  compreendido.  O  projeto  geralmente  desconsidera  efeitos  transientes  sobre 
o  desempenho  da  aeronave.  Tambem  se  costuma  assumir  uma  mudanga  de  forma  suave  entre  diferentes  configuragoes, 
o  que  permite  interpolagao  para  os  coeficientes  adimensionais  de  sustentagao  e  de  arrasto,  por  exemplo.  Todavia,  al¬ 
teragoes  podem  acontecer  num  curto  intervalo  de  tempo  em  avioes  de  combate,  por  exemplo.  Nestes  casos,  condigoes 
desfavoraveis  podem  resultar,  provocando  prejuizos  ao  sistema  de  controle  da  aeronave,  mesmo  que  momentaneos. 

Neste  sentido,  no  presente  trabalho  os  autores  investigam  o  comportamento  transiente  de  aerofolios  adaptativos.  Para 
tanto,  simulagoes  fluidodinamicas  computacionais  bidimensionais  sao  consideradas  para  calculo  dos  coeficientes  de  sus¬ 
tentagao,  de  arrasto  e  de  momento  de  arfagem  do  perfil  analisado.  Considera-se  um  aerofolio  NACA  0012  que  evolui  para 
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um  aerofolio  Selig  1223,  como  descrito  mais  a  frente. 

Para  facilitar  a  realizagao  das  simulagoes,  uma  metodologia  de  parametrizagao  espacial  de  aerofolios  foi  utilizada  para 
representagao  do  perfil  adaptativo;  uma  breve  descrigao  da  abordagem  considerada  e  apresentada  na  Se9ao  2.  Na  Segao  3 
sao  fornecidos  detalhes  das  simula9oes  numericas  conduzidas,  sao  apresentados  os  resultados  obtidos,  e  realizadas  as 
analises  pertinentes.  Por  fim,  conclusoes  sao  apresentadas  na  Se9ao  4. 

2.  PARAMETRIZAQAO  ESPACIAL  DE  AEROFOLIOS 

Para  realizagao  das  simulagoes  fluidodinamicas  transientes,  uma  unica  malha  estruturada  foi  concebida  para  a  situagao 
inicial  do  problema.  O  software  computacional  utilizado  (ANSYS  Fluent®,  versao  15.0)  foi  o  responsavel  por  promover 
as  alteragoes  geometricas  relacionadas  ao  aerofolio  adaptativo.  O  mesmo  dispoe  da  ferramenta  de  malha  dinamica,  que 
automaticamente  permite  suavizagao  (, smoothing )  ou  remalhagem  ( remeshing )  de  zonas  especificadas  pelo  usuario.  Para 
tanto,  uma  User  Defined  Function  (UDF)  -  um  codigo  personalizado  em  linguagem  C  -  deve  ser  implementada,  para  que 
o  software  realize  o  desejado  pelo  usuario.  Mais  a  respeito  deste  procedimento  sera  relatado  na  Segao  3. 

Com  respeito  a  UDF,  ela  e  a  responsavel  por  especificar  o  movimento  do  aerofolio  adaptativo  -  no  por  no.  Assim, 
um  esquema  de  parametriza9ao  espacial  foi  utilizado  para  representa9ao  do  aerofolio  adaptativo,  buscando  eliminar  com- 
plica9oes  relacionadas  a  altera9ao  de  malha  devidas  a  problemas  de  convergencia,  por  exemplo.  Desta  forma,  a  UDF 
implementada  torna-se  independente  das  pos^oes  ocupadas  pelos  nos  que  compoem  a  superficie  do  aerofolio  conside- 
rado.  Outra  vantagem  da  abordagem  considerada  e  permitir  a  otimiza9ao  da  forma  assumida  por  aerofolios  (Lane,  2010). 

A  estrategia  adotada  neste  trabalho  foi  proposta  por  Kulfan  e  Bussoletti  (2006).  A  representa9ao  matematica  utilizada 
e  a  seguinte: 


C  =  t"1  (1  -  *)Nl  Y  T,  +  t/sAg,  (1) 

i=o  lAn  l>' 

onde  £  =  z/c  Qif/  =  v/c,  sendo  c  a  corda.  O  par  (v,  z)  denota  um  ponto  do  aerofolio.  Ainda,  N\  e  N2  sao  constantes  que 
modificam  a  forma  da  fun9ao  dada  na  Eq.  (1),  e  Af  =  zte/c ,  sendo  zte  a  coordenada  do  aerofolio  em  seu  bordo  de  fuga, 
quando^  =  1.  Para  aerofolios  com  bordo  de  ataque  arrendondado  e  bordo  de  fuga  pontiagudo,  N\  =  0,5eA2  =  1.  Valores 
associados  a  outras  formas  geometricas  para  N\  e  N2  podem  ser  encontrados  na  literatura  (Kulfan,  2007,  2008).  Quanto 
aos  coeficientes  A/,  sao  os  responsaveis  por  dar  a  forma  particular  associada  a  um  aerofolio  especifico.  E  importante 
salientar  que  a  Eq.  (1)  deve  ser  aplicada  as  superficies  inferior  e  superior  de  um  aerofolio  de  forma  independente. 


Para  determina9ao  dos  valores  de  At ,  supoe-se  a  existencia  de  um  conjunto  de  pontos  (if/j  ,£/),  j  =  1 , . . . ,  N>n  + 1 ,  que 
descrevem  a  geometria  do  aerofolio  de  interesse.  Um  procedimento  de  ajuste  via  mmimos  quadrados  e  entao  considerado. 
Assim,  os  coeficientes  A;  constituem  a  solu9ao  do  sistema  linear  seguinte: 

(StS)A  =  ST(C  -  V>A£),  (2) 

onde: 

A  =  [A0  •••  Anf;  (3) 

Sj.M  =  [<'0  -  #j)N2]  ~  +jri ;  (4) 

C  =  [Ci  ■■■  £n]t ;  (5) 

ip  =  ■■■  i//n]T  ■  (6) 


Para  as  parametriza9oes  conduzidas  neste  trabalho,  tomou-se  n  =  20,  com  base  num  estudo  de  convergencia  do  erro 
medio  quadratico.  Ceze  et  al.  (2009)  alertam  que  valores  muito  grandes  de  n  tornam  o  sistema  linear  dado  na  Eq.  (2)  mal 
condicionado. 

Quanto  aos  pontos  (if/j,£j),  foram  obtidos  da  base  de  dados  da  University  of  Illinois  at  Urbana-Champaign  para  os 
aerofolios  de  interesse  (NACA  0012  e  Selig  1223).  Para  aumentar  o  numero  de  pontos  a  disposi9ao,  um  procedimento 
de  interpola9ao,  baseado  em  polinomios  de  Hermite  cubicos  (PCHIP  -  Piecewise  cubic  Hermite  interpolation ),  foi  tam- 
bem  considerado.  A  vantagem  do  metodo  de  parametriza9ao  adotado  com  respeito  a  este  procedimento  e  que  se  tern  a 
disposi9ao  expressoes  analiticas  para  as  curvas  que  compoem  um  dado  perfil,  validas  ao  longo  de  toda  sua  corda. 

3.  SIMULAQOES  NUMERICAS 

Nesta  Se9ao,  detalhes  sao  dados  sobre  o  cenario  assumido  para  simula9oes  numericas  e  sobre  os  metodos  utiliza- 
dos.  Embora  o  foco  principal  do  trabalho  envolva  as  caracteristicas  aerodinamicas  transientes  de  aerofolios  adaptativos, 
algumas  simula9oes  fluidodinamicas  em  regime  permanente  tambem  foram  conduzidas. 

O  software  ANSYS  Fluent®  versao  15.0  foi  aquele  escolhido  para  condu9ao  das  analises  apresentadas  a  seguir.  O 
modelo  de  turbulencia  empregado  foi  o  SST  (, Shear- Stress  Transport )  k-oj,  que  combina  a  formula9ao  padrao  do  modelo 
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k-co  para regioes proximas  a paredes,  onde  desempenha de  maneira robusta e  correta,  com  o  modelo  k-e  para regioes  de 
escoamento  livre  (Menter,  1994).  O  modelo  tambem  considera  efeitos  importantes  relacionados  ao  transporte  da  tensao 
de  cisalhamento  turbulenta  r  :=  - pu'v ',  o  que  melhora  seu  desempenho  quando  da  presenga  de  gradientes  adversos  de 
pressao  no  escoamento. 

O  valor  de  10-8  foi  assumido  para  verificagao  da  convergencia  dos  residuos  relacionados  as  equagoes  que  modelam 
o  escoamento.  O  metodo  de  solugao  acoplado  disponibilizado  pelo  software  ANSYS  Fluent®  foi  o  escolhido  para  lidar 
com  o  acoplamento  pressao-velocidade  das  equagoes  de  Navier- Stokes,  por  propiciar  melhores  taxas  de  convergencia.  Foi 
assumida  a  hipotese  de  incompressibilidade.  Gradientes  de  pressao  foram  interpolados  com  um  metodo  de  segunda  ordem. 
As  discretizagoes  para  os  gradientes  de  velocidade  e  das  variaveis  do  modelo  de  turbulencia  foram  computadas  com  o 
metodo  MUSCL  ( Monotone  Up  stream- Centered  Schemes  for  Conservation  Laws).  O  calculo  de  gradientes  no  centra 
de  celulas  foi  realizado  com  um  metodo  de  minimos  quadrados.  Quanto  ao  numero  maximo  de  iteragoes  consideradas 
para  convergencia  da  solugao,  este  foi  tornado  igual  a  250  por  passo  de  tempo  nas  simulagoes  transientes.  Em  regime 
permanente,  foram  realizados  calculos  ate  convergencia  dos  valores  observados  para  os  coeficientes  de  sustentagao,  de 
arrasto  e  de  momento  de  arfagem. 

Ainda,  as  simulagoes  foram  conduzidas  num  dommio  bidimensional.  A  corda  dos  aerofolios  considerados  foi  sempre 
tomada  igual  a  1  m,  e  o  angulo  de  ataque  mantido  com  o  valor  de  0°.  O  fluido  considerado  e  ar,  com  p  =  1,225  kg/m3  e 
p  =  1,7894  x  10-5  kg/(m-s).  Detalhes  sobre  condigoes  de  contorno  sao  fornecidos  na  Subsegao  3.2. 

3.1.  Fungoes  de  Interpolagao  Assumidas  para  Evolugao  da  Geometria  do  Aerofolio  Adaptativo 

Para  promover  a  alteragao  geometrica  do  aerofolio  adaptativo  quando  da  simulagao  fluidodinamica  transiente,  foram 
assumidas  evolugoes  de  tres  naturezas  distintas.  Essas  puderam  ser  especificadas  por  meio  de  uma  fungao  de  interpolagao 
i(t),  a  partir  da  qual  a  geometria  do  perfil  estudado  pode  ser  estabelecida  para  qualquer  instante  de  tempo,  por  meio  da 
equagao  seguinte: 

f  GM)  =  M)  i(0  +  fsOA)  [1  -  i(t)] ,  (7) 


onde  £n(*A)  e  £s(i//)  indicam  as  coordenadas  associadas  aos  perfis  NACA  0012  e  Selig  1223,  respectivamente. 

O  uso  da  Eq.  (7)  assume  que  a  corda  do  aerofolio  nao  se  altera  com  o  tempo.  Alem  do  mais,  algumas  condigoes  devem 
ser  respeitadas  pela  fungao  de  interpolagao,  sejam  elas:  (i)  que  i(t)  =  0  para  t  <  0  ;  e  (ii)  que  i(t)  -  1  para  t  >  tf ,  onde  ti 
etf  representam  os  instantes  de  inicio  e  de  fim  de  alteragao  geometrica  do  aerofolio  adaptativo. 

As  expressoes  das  fungoes  de  interpolagao  assumidas  neste  estudo,  para  ti  <  t  <  tf,  sao  fornecidas  a  seguir: 


tf  ti 


t-L  1  /  t-U 

^harmonica  —  v 


^'polinomial  =  35 


t-tj 
tf  ~  ti 


-84 


tf  ti 

t-tj 
tf  -  ti 


+  70 


t-tj 
tf  -  ti 


-20 


t-tj 
tf  -  ti 


(9) 

(10) 


A  Fig.  1  ilustra  o  comportamento  de  cada  uma  das  anteriores  no  intervalo  ti  <  t  <  tf,  enquanto  a  Fig.  2  mostra  a 
evolugao  do  contorno  do  aerofolio  para  o  caso  em  que  interpolagao  linear  e  considerada.  Cabe  destacar  que,  enquanto 
a  fungao  de  interpolagao  linear  apresenta  derivada  primeira  descontmua  para  t  =  ti  e  t  =  tf,  as  fungoes  harmonicas  e 
polinomial  sao  mais  suaves,  utilizadas  inclusive  para  o  desenvolvimento  de  sistemas  de  cames  (Waldron  e  Kinzel,  1999). 
A  fungao  de  interpolagao  harmonica  apresenta  derivada  de  terceira  ordem  com  descontinuidade  em  t  =  ti  e  t  =  tf,  e 
aquela  polinomial  tern  derivada  de  ordem  quatro  descontmua. 

As  diferengas  quanto  a  continuidade  de  derivadas  e  o  que  motivou  a  adogao  das  diferentes  fungoes  de  interpolagao 
apresentadas  anteriormente.  A  aplicagao  brusca  de  esforgos  pode  acarretar  vibragoes  e  oscilagoes  transientes  na  superficie 


^  I  inear 
l  harmonica 
1  polinomial 


Figura  1.  Comportamento  das  fungoes  de  interpolagao  consideradas. 
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do  aerofolio  num  caso  pratico.  Assim,  mesmo  apos  o  periodo  de  alteragao  geometrica  do  aerofolio,  podem  resultar  flutu- 
agoes  nos  coeficientes  de  sustentagao,  de  arrasto  e  de  momento  de  arfagem.  Como  consequencia,  um  sistema  de  atuagao 
mais  robusto  seria  desejado,  com  vistas  a  suplantar  eventuais  complicagoes  introduzidas  durante  etapa  de  modificagao  de 
forma.  Como  se  comprovara  com  a  apresentagao  dos  resultados,  pode  ser  desejavel  a  adogao  de  uma  estrategia  especifica 
para  aplicagao  dos  esforgos  de  atuagao  a  estrutura  de  um  dado  aerofolio  adaptativo,  para  melhora  de  suas  caracteristicas 
transientes. 

3.2.  Geragao  de  Malhas  e  CondigSes  de  Contorno 

As  malhas  consideradas  neste  estudo  foram  obtidas  com  uso  do  software  ANSYS ICEM  CFD®  versao  15.0.  Devido  as 
analises  de  regime  permanente  conduzidas,  nove  malhas  distintas  foram  produzidas.  Todas  contam  com  aproximadamente 
o  mesmo  numero  de  celulas  e  de  nos,  por  volta  de  200  000,  e  apresentam  qualidade  superior  a  90%.  A  regiao  proxima  ao 
aerofolio  das  malhas  dos  aerofolios  NACA  0012  e  Selig  1223  podem  ser  visualizadas  na  Fig.  3. 

Independentemente  do  aerofolio  considerado,  a  estrategia  de  blocagem  adotada  foi  do  tipo  C.  A  distancia  assumida 
para  aplicagao  de  condigoes  de  contorno  foi  de  30  cordas.  Para  a  primeira  camada  de  elementos  junto  a  um  dado  aerofolio 
foi  considerada  uma  espessura  de  10-6  m,  o  que  garante  modelagem  de  camadas  limites  para  as  analises  conduzidas 
(y+<0,01  para  a  condigao  de  entrada  especificada).  Os  elementos  proximos  as  fronteiras  externas  do  dominio  tern  dimen- 
sao  caracteristica  aproximadamente  igual  a  0,2  m. 

Quanto  as  condigoes  de  contorno,  foram  consideradas:  (i)  de  nao  deslizamento  e  de  nao  penetragao  junto  ao  aerofolio; 
(ii)  de  simetria  nas  faces  superior  e  inferior  da  malha;  (iii)  de  saida  na  face  a  jusante,  com  pressao  nula;  e  (iv)  de  entrada, 
com  velocidade  igual  a  1,460735  m/s  ao  longo  da  diregao  horizontal  na  face  a  montante  do  aerofolio.  Resulta  um  numero 
de  Reynolds  igual  a  105.  Tambem  resulta  que  o  numero  de  Mach  tern  valor  baixo,  suficiente  para  justificar  a  hipotese  de 
incompressibilidade.  Com  respeito  a  turbulencia,  especificou-se  intensidade  turbulenta  (definida  como  a  razao  do  valor 
RMS  de  flutuagoes  com  respeito  a  velocidade  media  do  escoamento)  igual  a  5%,  e  razao  de  viscosidade  turbulenta  (nt  In) 
igual  a  10  tanto  na  entrada  quanto  na  saida. 

3.3.  Caracteristicas  Aerodinamicas  dos  Aerofolios  Interpolados  em  Regime  Permanente 

Simulagoes  de  regime  permanente  foram  realizadas  para  obtengao  das  caracteristicas  operacionais  relacionadas  a 
alguns  dos  aerofolios  interpolados,  observados  durante  a  transigao  do  aerofolio  adaptativo  considerado  neste  estudo. 

A  partir  de  tais  simulagoes,  os  valores  do  coeficiente  de  sustentagao,  Q ,  do  coeficiente  de  arrasto,  Cj,  e  do  coeficiente 
de  momento  de  arfagem  medido  a  1/4  do  valor  da  corda  a  partir  do  bordo  de  ataque,  (Cm)i/4,  puderam  ser  estimados. 
Para  efeito  de  comparagao,  analises  tambem  foram  feitas  com  uso  do  software  XFOIF®  versao  6.99.  Este  codigo  faz 
uso  de  um  metodo  de  paineis  linear  com  corregoes  para  efeitos  oriundos  de  compressibilidade  para  calculos  nao  viscosos 


a)  NACA  0012 


b)  Selig  1223 


Figura  3.  Detalhe  de  duas  das  malhas  construidas,  proximas  ao  aerofolio. 
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(Drela,  1989).  Fontes  distribuidas  sao  superpostas  ao  longo  do  aerofolio  e  da  esteira  para  compensar  efeitos  da  ultima 
natureza.  Equagoes  integrals  para  camadas  limite  sao  resol vidas  em  conjunto  com  as  equagoes  que  modelam  o  escoamento 
inviscido.  Apesar  de  ser  uma  ferramenta  que  nao  propicia  resultados  com  alta  fidelidade,  e  rapida,  e  por  este  motivo  e 
comumente  empregada  em  etapas  preliminares  de  analises  de  compromisso. 

Os  graficos  mostrados  na  Fig.  4  evidenciam  que  Q  e  (Cm)  1/4  alteram-se  de  acordo  com  a  lei  de  interpolagao  utilizada 
na  analise,  /linear-  Quanto  ao  comportamento  de  Cd,  este  e  distinto.  Credita-se  estas  observagoes  ao  fato  de  Q  e  de 
(Cm)  1/4  dependerem principalmente  da  distribuigao  de  pressao,  e  de  Cd  ser  fungao  das  tensoes  cisalhantes  junto  a  parede, 
para  o  angulo  de  ataque  considerado.  Alem  disso,  com  o  aumento  de  arqueamento,  a  medida  que  o  aerofolio  adaptativo 
aproxima-se  daquele  Selig  1223,  uma  maior  regiao  de  separagao  encontra-se  presente  junto  ao  bordo  de  fuga,  o  que 
explica  o  aumento  de  Cd  da  forma  como  ocorre. 

Observam-se  ainda  distingoes  entre  os  resultados  obtidos  atraves  do  softwares  ANSYS  Fluent®  e  XFOIL®,  em  espe¬ 
cial  para  Cj.  Estas  diferengas  sao  creditadas  ao  baixo  numero  de  Reynolds  considerado,  e  as  zonas  de  recirculagao  no 
extradorso  e  intradorso  do  aerofolio  Selig  1223,  que  tendem  a  ocorrer  para  aerofolios  interpolados  quando  /linear  >0,5. 
Cabe  destacar  que  os  valores  estimados  pelo  software  ANSYS  Fluent®  para  C/,  Cd  e  (Cm)  1/4  encontram-se  proximos 
daqueles  observados  em  experimentos  conduzidos  com  os  perfis  NACA  0012  (Sheldahl  e  Klimas,  1981)  e  Selig  1223 
(Selig  et  al. ,  1995),  ainda  que  as  simulagoes  realizadas  tenham  sido  bidimensionais. 

3.4.  User  Defined  Functions 

Para  cumprir  com  a  alteragao  geometrica  do  aerofolio  adaptativo  considerado,  uma  User  Defined  Function  foi  desen- 
volvida.  Como  ja  mencionado,  uma  UDF  trata-se  de  um  codigo  computacional  em  linguagem  C.  O  software  ANSYS 
Fluent®  dispoe  de  tal  funcionalidade  para  atender  diversas  necessidades  do  usuario,  que  podem  envolver  a  especificagao 
de  uma  condigao  de  contorno  nao  uniforme,  ou  a  promogao  de  alteragoes  geometricas,  por  exemplo. 

O  codigo  desenvolvido  em  linguagem  C  conta  com  macros  que  facilitam  a  implementagao  de  uma  UDF.  Para  o  caso 
em  questao,  o  software  ANSYS  Fluent®  ja  disponibiliza  ao  usuario  macros  utilizadas  para  fazer  lagos  sobre  nos  de  uma 
dada  zona  da  malha,  assim  como  para  especificar  a  posigao  de  nos  de  forma  individual.  Desta  forma,  o  codigo  desenvol¬ 
vido  neste  trabalho  apenas  implementa  as  alteragoes  relacionadas  ao  aerofolio  adaptativo  considerado,  com  auxilio  das 
equagoes  mostradas  anteriormente.  Todas  as  operagoes  realizadas  sobre  a  malha  sao,  a  partir  dai,  feitas  com  auxilio  de 
macros  disponibilizadas  pelo  software. 

Quanto  a  implementagao  da  UDF  junto  ao  software  ANSYS  Fluent®,  para  o  caso  em  questao,  esta  deve  ser  feita 
por  meio  de  compilagao  do  codigo  no  ambiente  do  proprio  software.  Cabe  destacar  que  tal  compilagao  e  dependente  da 
plataforma  operacional  utilizada,  e  esta  sujeita  a  disponibilidade  de  fungoes  de  suporte  no  caminho  do  sistema,  as  quais 
podem  ser  encontradas  em  pacotes  de  desenvolvedores  de  software.  Depois  de  compilada,  a  UDF  deve  ser  carregada  no 
ambiente  do  software  ANSYS  Fluent®,  a  partir  do  que  se  encontra  dispomvel  para  uso  por  meio  da  interface  grafica  do 
programa  atraves  de  menus  apropriados. 

Com  respeito  aos  desenvolvimentos  feitos  neste  trabalho,  pode-se  mencionar  que  o  principal  empecilho  para  a  utili- 
zagao  de  UDFs  que  promo vem  alteragao  geometrica  de  malha  esta  relacionado  a  criagao  de  celulas  com  areas/volumes 
negativos  entre  um  passo  de  tempo  e  outro  da  simulagao.  Muitas  das  vezes  nao  se  pode  identificar  onde  ocorre  o  pro- 
blema,  e  a  redugao  do  passo  de  tempo,  ou  o  estabelecimento  de  criterios  mais  rigorosos  para  suavizagao  da  malha  ou 
remalhagem,  podem  nao  ser  suficientes  para  contornar  a  situagao. 

Outro  aspecto  relacionado  a  alteragao  da  geometria  do  aerofolio  adaptativo  diz  respeito  a  modificagao  da  malha. 
Aspectos  teoricos  relacionados  aos  procedimentos  adotados  pelo  software  ANSYS  Fluent®  podem  ser  encontrados  em 
sua  documentagao  (ANSYS,  2013).  A  metodologia  preve  remalhagem  ou  suavizagao  da  malha  por  meio  de  movimentagao 
de  seus  nos.  Para  os  casos  apreciados  neste  trabalho,  a  segunda  estrategia  foi  considerada  por  estar  associada  a  menores 
custos  computacionais.  A  qualidade  da  malha  quando  de  sua  suavizagao  pode  ser  garantida  junto  a  parede  do  aerofolio 


Figura  4.  Caractensticas  aerodinamicas  dos  aerofolios  interpolados  para  a  condigao  de  regime  permanente. 
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pela  especificagao  de  uma  regiao  de  camada  limite,  que  evita  distorgoes  nas  celulas,  e  a  mantem  ortogonal  a  superficie 
durante  todo  o  tempo  de  simulagao  computacional,  o  que  contribui  para  a  obtengao  de  resultados  mais  fidedignos. 

3.5.  Resultados  Transientes 

Os  resultados  transientes  apresentados  a  seguir  foram  obtidos  com  o  auxilio  de  uma  UDF  e  com  condigoes  iniciais  e 
de  contorno  ja  detalhadas  ao  longo  do  texto. 

Para  todos  os  casos  considerados,  tomaram-se  os  tempos  de  imcio  e  de  fim  de  alteragao  geometrica  do  perfil  adaptativo 
iguais  a  1  s  e  2  s,  respectivamente.  O  tempo  final  para  cada  uma  das  simulagoes  foi  escolhido  igual  a  10  s,  para  permitir 
observagao  da  evolugao  das  caracteristicas  operacionais  do  aerofolio  depois  de  passado  o  periodo  de  atuagao.  O  passo  de 
tempo  para  integragao  foi  mantido  fixo,  igual  a  0,001  s.  Um  metodo  implicito  de  segunda  ordem  foi  considerado. 

As  Figs.  5,  6  e  7  mostram  as  evolugoes  de  Q,  Cd  e  (Cm)  1/4,  respectivamente,  para  o  intervalo  de  tempo  considerado. 
Pode-se  constatar  por  meio  destas  que  uma  alteragao  de  forma  linear  esta  relacionada  a  mudangas  abruptas  nos  valores 
dos  coeficientes  adimensionais  considerados,  durante  o  imcio  e  o  fim  da  reconfiguragao  do  aerofolio.  Alem  disso,  ha 
um  salto  nos  valores  dos  coeficientes  nestes  instantes  especificos  de  tempo.  Quando  a  interpolagao  e  harmonica,  uma 
transigao  mais  suave  e  obtida,  de  forma  que  saltos  nao  estao  presentes  nos  historicos  apresentados.  Contudo,  como  se 
observa  nas  Figs.  5c  e  5d,  as  derivadas  de  Q  sao  descontmuas  durante  o  imcio  e  o  fim  da  transigao  do  perfil  analisado, 
o  que  tambem  ocorre  para  Cj  e  (Cm)  1/4.  Uma  das  razoes  para  o  surgimento  de  tais  descontinuidades  diz  respeito  aos 
algoritmos  utilizados  pelo  software  ANSYS  Fluent®.  Os  autores,  contudo,  asseguram  a  convergencia dos  residuos  a  mveis 
apropriados,  de  forma  a  permitir  a  analise  dos  resultados  apresentados.  Nos  casos  de  interpolagao  linear  ou  harmonica,  as 
descontinuidades  associadas  as  derivadas  dos  coeficientes  aerodinamicos  considerados  levariam  a  aplicagao  de  esforgos 
tipo  impacto  sobre  a  estrutura  do  aerofolio.  Tais  forgas  ocasionariam  potenciais  problemas  do  ponto  de  vista  de  controle 


s  T  ,  •  A  t  .  ~  Tempo  [s] 

c)  Imcio  da  transigao 


Figura  5.  Historico  do  coeficiente  de  sustentagao  relacionado  ao  aerofolio  adaptativo  considerado. 


Figura  6.  Historico  do  coeficiente  de  arrasto  relacionado  ao 
aerofolio  adaptativo  considerado.  (Legenda  como  na  Fig.  5.) 
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b)  Transigao 


Tempo  [s] 


Figura  7.  Historico  do  coeficiente  de  momento  de  arfagem  relacionado 
ao  aerofolio  adaptativo  considerado.  (Legenda  como  na  Fig.  5.) 

da  aeronave  durante  uma  etapa  seguramente  critica,  de  reconfiguragao  geometrica.  Por  este  motivo,  o  sistema  de  atuagao 
que  confere  adaptatividade  ao  aerofolio  necessitaria  de  elevada  robustez,  mesmo  que  apenas  para  nao  sofrer  danos,  e 
evitar  sua  posterior  falha,  frente  a  ocorrencia  das  descontinuidades  constatadas  nos  esforgos  aerodinamicos. 

A  implementagao  de  uma  alteragao  geometrica  por  meio  de  interpolagao  polinomial,  por  outro  lado,  favorece  uma 
transigao  adequada,  no  sentido  de  exibir  initio  e  fim  caracterizados  por  um  comportamento  suave.  Apesar  deste  fato, 
nenhuma  das  interpolates  consideradas  para  promover  a  evolugao  transiente  do  aerofolio  adaptativo  considerado  permite 
a  obtengao  de  caracteristicas  aerodinamicas  proximas  aquelas  observadas  em  regime  permanente.  Alem  disso,  mesmo 
passados  8  s  depois  do  fim  da  reconfiguragao  geometrica  do  aerofolio,  regime  permanente  ainda  nao  foi  alcangado,  como 
comprovado  pelas  Figs.  5,  6  e  7  -  o  valor  de  Q,  por  exemplo,  tende  a  aumentar  ate  um  patamar  pouco  superior  a  1,0. 

Para  complementar  os  resultados  apresentados,  as  Figs.  8,  9  e  10  mostram  os  vetores  velocidades,  coloridos  de  acordo 
com  sua  magnitude,  proximos  a  superficie  do  aerofolio  adaptativo,  para  os  casos  de  interpolagao  linear,  harmonica  e 
polinomial,  respectivamente.  Pode-se  observar  a  formagao  de  uma  bolha  no  intradorso  do  aerofolio  para  os  tres  casos 
considerados.  A  mesma  se  propaga  pela  referida  superficie  apos  o  instante  em  que  tern  fim  a  etapa  de  alteragao  geometrica. 
Quando  do  estendimento  de  tal  bolha  ate  o  bordo  de  fuga  do  aerofolio,  ocorre  um  decrescimo  momentaneo  nos  coeficientes 
de  sustentagao  e  de  momento  de  arfagem,  e  uma  grande  queda  no  valor  do  coeficiente  de  arrasto,  quando  o  tempo  vale 
aproximadamente  4  s.  Um  pequeno  atraso  observado  na  evolugao  de  Q  para  o  caso  de  interpolagao  linear,  como  mostra  a 
Fig.  5a,  pode  ser  atribuido  a  formagao  de  uma  bolha  de  menor  tamanho  ao  final  da  etapa  de  transi^o  geometrica,  quando 


Figura  8.  Vetores  velocidade  relacionados  a  evolu^ao  do  perfil  adaptativo  segundo  interpola^ao  linear. 
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Figura  9.  Vetores  velocidade  relacionados  a  evolu^ao  do  perfil  adaptativo  segundo  interpola^ao  harmonica. 


/  =  1,4  s 


Figura  10.  Vetores  velocidade  relacionados  a  evolu^ao  do  perfil  adaptativo  segundo  interpola^ao  polinomial. 
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sao  comparadas  as  figuras  que  ilustram  os  vetores  velocidade  relacionados  as  simulagoes  realizadas.  Outro  ponto  que 
merece  destaque  e  a  semelhanga  entre  os  resultados  relacionados  as  condigoes  de  interpolagao  harmonica  e  polinomial. 
As  poucas  diferengas  observadas  se  refletem  num  comportamento  transiente  muito  proximo  dos  coeficientes  Q ,  Cd  e 

(Cm)l/4- 

4.  CONCLUSOES 

O  presente  trabalho  buscou  esclarecer  aspectos  relacionados  ao  comportamento  transiente  exibido  por  aerofolios  adap- 
tativos  atraves  de  simulagoes  computacionais  bidimensionais  de  alta  fidelidade.  Neste  contexto,  um  codigo  em  linguagem 
C  -  UDF  -  foi  desenvolvido  com  o  intuito  de  introduzir  as  alteragoes  geometricas  cabiveis  no  ambiente  do  software  con- 
siderado  para  simulagoes.  Uma  metodologia  de  parametrizagao  espacial  de  aerofolios  foi  ainda  considerada  para  facilitar 
a  implementagao  do  referido  codigo. 

Algumas  criticas  podem  ser  feitas  a  metodologia  considerada.  A  principio,  em  uma  situagao  pratica,  o  sistema  de 
atuagao  de  um  aerofolio  adaptativo  muito  dificilmente  conseguiria  promover  as  alteragoes  geometricas  imaginadas  neste 
estudo,  seja  pela  grande  distingao  das  formas  associadas  aos  estados  inicial  e  final  do  aerofolio,  ou  pela  energia  de  atuagao 
requerida  para  promover  tal  mudanga.  Entretanto,  em  casos  praticos,  a  metologia  pode  ser  aplicada  apos  a  identificagao  do 
historico  real  de  evolugao  da  geometria  do  aerofolio,  considerando  as  limitagoes  dos  atuadores.  Outro  ponto  que  contribui 
para  a  nao  execugao  de  uma  transigao  como  aquelas  idealizadas  diz  respeito  a  interagao  fluido-estrutura,  que  pode  ajudar  a 
evitar  descontinuidades  percebidas  nos  resultados  computados,  ou  entao  levar  a  problemas  de  natureza  aeroelastica,  como 
divergencia  e  flutter. 

Independentemente  destas  restrigoes,  os  resultados  apresentados  sao  conclusivos  no  sentido  de  estabelecer  uma  relagao 
entre  a  forma  de  aplicagao  de  esforgos  a  estrutura  do  aerofolio,  e  seu  comportamento  aerodinamico  durante  o  transiente  de 
sua  evolugao.  Neste  sentido,  pode-se  estabelecer  como  objetivo  de  uma  eventual  otimizagao  o  cumprimento  de  metricas 
de  desempenho  durante  o  regime  de  reconfiguragao  do  aerofolio  adaptativo.  Desta  forma,  pode-se  evitar  a  ocorrencia 
de  descontinuidades,  e  promover  uma  evolugao  de  C/,  Cd  e  (Cm)  1/4  que  evite  sobrecarregar  o  sistema  de  controle  da 
aeronave,  por  exemplo. 

Tambem  fica  evidente  que  a  evolugao  transiente  de  aerofolios  adaptativos  nao  valida  a  pratica  de  interpolagao  de 
resultados  obtidos  atraves  de  analises  que  consideram  regime  permanente.  Efeitos  de  natureza  dinamica  sao  importantes,  e 
ocasionam  um  atraso  significativo  ate  que  eventual  regime  estacionario  seja  alcangado.  Assim,  o  desempenho  do  aerofolio 
adaptativo  apresenta  um  comportamento  particular  durante  sua  transigao,  que  e  dependente  da  forma  de  atuagao  -  e 
consequente  estrategia  de  mudanga  geometrica. 

Varios  aspectos  relacionados  a  este  trabalho  ainda  necessitam  investigagao,  e  devem  ser  considerados  em  estudos 
futuros.  Entre  eles  se  encontram:  (i)  a  influencia  do  tempo  total  gasto  para  transigao  geometrica;  (ii)  a  influencia  do 
numero  de  Reynolds  sobre  as  caracteristicas  transientes  do  aerofolio  adaptativo;  e  (iii)  a  consideragao  de  efeitos  de  com- 
pressibilidade  e  tridimensionais  sobre  o  comportamento  aerodinamico  do  aerofolio  adaptativo.  Ainda  e  de  interesse  dos 
autores  avaliar  estrategias  alternativas  para  previsao  de  escoamentos  sobre  aerofolios  com  geometria  variavel  (Gaunaa, 
2006;  Walker,  2009;  Johnston  et  al. ,  2010). 
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Abstract:  Morphing  aircraft  have  been  extensively  studied  due  to  their  promised  benefits.  For  instance,  it  is  claimed 
that  morphing  aircraft  are  able  to  diminish  fuel  consumption  due  to  their  greater  aerodynamic  efficiency  throughout  the 
aircraft  flight  envelope.  Although  several  research  studies  have  investigated  morphing  aircraft  design,  which  have  been 
followed  by  prototype  manufacturing  and  tests,  few  have  focused  on  the  transient  behavior  that  results  from  the  morphing 
process  itself.  In  this  paper,  one  presents  the  aerodynamic  characteristics  (lift,  drag  and  moment  coefficients)  related 
to  an  adaptive  airfoil,  that  initially  assumes  the  geometry  of  a  NACA  0012  profile,  and  evolves  to  a  Selig  1223  profile. 
Two-dimensional  simulations  of  the  flow  surrounding  the  airfoil  are  performed  with  the  ANSYS  Fluent®  15.0  software.  A 
spatial  parameterization  scheme  is  adopted  to  implement  the  necessary  geometric  shape  changes.  Results  are  compared 
with  steady  state  simulations’  data.  They  show  that  the  history  of  the  geometric  shape  change  has  a  rather  great  influence 
on  the  computed  aerodynamic  coefficients,  thus  leading  to  eventual  discontinuities  and  issues  for  the  aircraft  control 
subsystem. 

Keywords:  Morphing  Aircraft,  Computational  Fluid  Dynamics. 
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Abstract.  This  paper  aims  at  the  investigation  of  the  behavior  of  the  flow  over  a  low  frequency  pitching  NACA  0012 
airfoil  at  Reynolds  number  Re  =  100,000  by  the  analysis  of  a  two-dimensional  RANS  simulation.  The  behavior  of  the 
flow  wake  at  the  trailing  edge  is  studied  by  the  analysis  of  the  vorticity  for  each  incidence  angle  and  results  are 
compared  by  the  study  of  theoretical  concepts  and  experimental  data.  The  use  of  standard  Reynolds -Averaged  Navier- 
Stokes  (RANS)  simulation  has  shown  accuracy  on  the  prediction  of  stall  and  reattachment  incidence  angles  for 
upstroke  and  downstroke.  The  simulations  were  also  capable  of  capturing  flow  information  that  is  in  agreement  with 
experiments  despite  the  overprediction  of  lift  and  drag  coefficients,  due  to  the  two-dimensional  approximation  used  for 
the  analysis.  The  study  has  also  compared  the  two-dimensional  RANS  SST  model  simulation  data  with  previous  results 
from  three-dimensional  simulation  Delayed  Detached  Eddy  Simulation  (DDES)  SST  model  analysis,  showing  that  in 
comparison  of  the  three  components  analysis  the  approximation  captures  accurately  same  flow  characteristics.  The 
simulations  also  confirm  the  Coanda  Effect,  where  the  influence  of  the  lower  surface  of  the  airfoil  has  a  strong 
influence  due  to  flow  attachment  and  bending  in  the  trailing  edge,  acting  together  with  the  vortices  in  the  upper  surface 
in  the  lift  loss  of  the  dynamic  stall. 

Keywords:  Pitching  Airfoil,  Flapping  Wings,  Unsteady  Flow,  Low  Reynolds  Number,  Low  Oscillatory  Frequency. 


1.  INTRODUCTION 

The  flow  past  airfoils  has  been  studied  for  almost  a  century  when  considering  the  Thin  Airfoil  Theory  (Munk, 
1922)  applied  to  steady  airfoils.  These  studies  are  the  edge  of  aeronautical  research,  which  allows  the  ongoing 
optimization  of  aeronautical  profiles  in  modern  aircraft.  Unsteady  effects  used  to  be  ignored  for  simplicity  during 
experimental  research;  hence  the  flow  over  conventional  fixed  airfoils  is  widely  studied  and  fairly  well  understood. 
However,  this  approximation  is  not  sufficient  to  model  the  turbulent  flow  in  the  trailing  edge  during  flight.  Unsteadiness 
is  an  inherent  part  of  the  flow  during  flight,  where  there  are  far  more  variables  that  can  change  and  cannot  be  accounted 
for  during  wind  tunnel  steady  studies  made  on  airfoils.  These  studies  used  to  consider  a  well -developed  flow  past  the 
airfoils  to  predict  static  stall  and  aerodynamics  coefficients  as  variables  of  the  angle  of  incidence. 

For  steady  airfoils,  the  flow  encounters  the  airfoil  wall  and  attaches  to  it,  forming  a  laminar  boundary  layer.  In  the 
transitioning  from  laminar  flow  to  turbulent  flow,  a  separation  bubble  is  formed  close  to  the  airfoil  leading  edge, 
creating  a  reverse  flow  region.  Afterwards  the  boundary  layer  is  reattached  into  turbulent  flow  until  it  reaches  the 
trailing  edge  turbulent  separation  point,  and  thus  creates  a  separated  turbulent  shear  layer  and  the  detached  turbulent 
separation  region.  As  the  incident  angles  of  the  airfoil  increases,  the  trailing  edge  separation  point  progress  upward  in 
the  upper  surface  until  it  reaches  the  transition  bubble.  At  this  point  the  flow  does  not  reattach  after  the  laminar 
separation,  the  bubble  “bursts”  and  then  a  separated  turbulent  flow  region  is  created,  leading  to  the  stall  of  the  airfoil  by 
losing  lift.  The  flow  behaves  in  a  different  way  for  unsteady  airfoils,  since  the  vertical  wake  is  now  time -dependent  due 
to  the  unsteadiness  and  the  aerodynamics  coefficients  can  change  accordingly.  The  stall  of  a  body  under  unsteady 
motion  is  quite  complex  when  compared  to  static  stall  (McCroskey,  1982). 

The  studies  on  unsteady  airfoils  first  began  in  the  end  of  the  1970’s,  as  an  attempt  to  model  and  understand  the 
behavior  of  the  flow  when  considering  the  addition  of  the  fluctuation  of  both  body  and  flow  and  to  study  the  turbulence 
in  the  trailing  edge  flow.  Due  to  difficulties  on  data  acquisition  in  the  experiments  and  lack  of  computing  capability  for 
computational  turbulence  modelling,  the  studies  on  the  field  were  stopped  for  the  past  fifteen  years.  Advancements 
made  on  data  acquisition  devices  in  the  end  of  the  1990’s,  as  PIV  (Particle  Image  Velocimetry)  and  LDV/LDA  (Laser 
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Doppler  Velocimeter/ Anemometer),  were  essential  to  improve  experimental  measurements  and  wake  studies,  as  shown 
by  Berton  et  al.  (1997).  Also,  the  advancements  on  computing  memory  operations  and  development  of  URANS 
methods,  as  the  k-co  SST  model  (Menter,  1993),  and  LES  methods,  were  essential  to  the  turbulence  modelling  of 
unsteady  airfoils. 

Most  aeronautical  devices  can  encounter  this  situation  during  flight,  whether  designed  for  that  or  not.  The  unsteady 
flow  is  easily  observable  in  the  flow  of  rotorcraft  devices  and  wind  turbines,  but  also  it  is  present  in  fixed  wings,  which 
can  vibrate  at  high  frequencies  during  flight,  and  in  the  most  fascinating  spectacle,  the  flight  of  birds,  which  are  not  yet 
well  understood  due  to  the  challenges  in  creating  and  simulating  a  device  that  can  reproduce  all  the  motion 
characteristics  of  bird’s  flapping  wings.  Nowadays,  the  modelling  of  flapping  wings,  or  flapping  airfoils  when 
considering  a  two  dimensional  approach,  is  the  combination  of  a  plunging  (or  heaving)  and  pitching  (oscillatory) 
motion.  When  considering  this  kind  of  motion,  it  is  possible  to  see  beneficial  effects  of  unsteadiness,  which  is 
substantially  important  to  the  propulsive  efficiency  of  flapping  motion.  The  unsteadiness  is  also  a  way  to  delay  dynamic 
stall  to  control  periodic  vortex  generation  and  improve  the  performance  of  rotorcrafts  and  wind  turbines  (McCroskey, 
1982). 

As  a  component  of  a  flapping  motion,  a  pitching  airfoil  is  a  simple  way  to  study  and  can  be  used  to  assess  the 
influence  of  the  trailing  edge  vortex  on  flow  reattachment.  Due  to  the  flow  fluctuation,  the  circulation  near  the  wall 
varies,  which  is  accompanied  by  the  shedding  of  free  vorticity  into  the  wake.  When  the  pitching  airfoil  reaches  a  high 
incidence  angle,  the  vortex  has  a  high  energy  profile,  causing  vortex  shedding  at  the  leading  edge  and  reaching  dynamic 
stall,  causing  large  loss  in  the  lift  coefficient  and  increasing  drag.  As  the  incidence  angle  decreases  downstroke,  the  flow 
reattaches,  delayed  as  compared  to  static  stall,  as  found  in  the  experiments  of  McAlister  et  al.  (1978)  and  Berton  et  al. 
(2002). 

2.  TURBULENCE  MODELLING 


Aeronautical  flows  are  surely  a  class  where  the  prediction  of  its  properties  needs  high  accuracy,  mainly  due  to 
strong  adverse  pressure  gradients  and  separation  in  boundary  layers.  The  Reynolds -averaged  Navier-Stokes  equations 
(or  RANS  equations)  are  time -averaged  equations  of  motion  for  fluid  flow  and  can  be  used  with  approximations  based 
on  knowledge  of  the  properties  of  flow  turbulence  to  give  approximate  time-averaged  solutions  to  the  Navier-Stokes 
equations.  The  flow  solver  is  the  Fluent  14.5  (ANSYS,  Inc.,  2012). 

The  k-s  and  k-co  two-equation  RANS  models  are  not  able  to  capture  the  proper  behavior  of  turbulence  in 
aeronautical  flow.  The  popular  k-s  can  give  a  well-defined  boundary  layer-edge  during  simulation;  however,  it  is  less 
accurate  and  complex  on  sublayer  modelling.  The  k-co  is  substantially  more  accurate  in  the  sublayer;  yet  it  is  sensitive 
in  the  freestream,  which  is  the  cause  of  the  k-s  being  the  standard  equation  in  turbulence  modelling  (Menter  et  al ., 
2013).  Both  standard  two-equation  models  overpredict  the  shear  stress  in  adverse  pressure  gradient  flows,  even  when 
considering  delayed  separation.  The  Shear  Stress  Transport  SST  model  (Menter,  1993)  was  developed  due  to  the  need 
of  more  accurate  separation  prediction  for  aeronautic  flows.  The  k-co  SST  model  is  a  blend  of  a  k-co  model,  which  is 
used  near  walls  in  the  sublayer  prediction,  and  a  k-s  model,  used  to  predict  the  flow  in  the  freestream  region.  Thus,  the 
model  is  fairly  robust,  since  it  can  accurately  predict  the  flow  at  both  sublayer  and  boundary  layer  edge  and  works  better 
at  capturing  recirculation  regions  by  enforcing  the  Bradshaw  Relation. 

To  blend  the  k-s  and  k-co  models,  it  is  necessary  to  transform  the  former  into  equations  based  on  k  and  co.  This  leads 
to  the  cross-diffusion  term,  defined  in  Eq.  1 . 
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The  k-co  SST  Turbulence  Kinect  Energy  function  is  given  as  shown  in  Eq.  (2).  The  Dissipation  Rate,  combining 
both  standard  k-co  e  k-s  models  by  the  use  of  the  cross-diffusion  term  and  the  blending  function  is  shown  at  Eq.  (3). 
There  are  the  two  equations  of  this  class  of  RANS  model. 
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The  blending  function  is  defined  by  Eq.  (4),  which  is  equal  to  zero  in  the  freestream  region,  activating  the  k-s  cross 
diffusion  term,  and  switches  over  to  one  in  the  boundary  layer  zone,  to  assure  an  accurate  calculation  by  the  use  of  the 
k-co  function. 
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The  Turbulent  Eddy  Viscosity  for  the  k-co  SST  model  is  defined  as  follows: 


(4) 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


VIM  Congresso  Nacional  de  Engenharia  Mecanica,  10  a  15  de  agosto  de  2014,  Uberlandia  -  Minas  Gerais 
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(5) 


This  term  enforces  the  Bradshaw  relation,  to  guarantee  accurate  separation  prediction  for  the  model.  The  Bradshaw 
relation  is  given  as: 

— uv  =  vt  |^|  =  axk  (6) 


Back  to  the  Turbulent  Eddy  Viscosity  equation,  S  is  the  invariant  measure  of  the  strain  rate,  while  F2  is  a  second 
blending  function,  as  shown  in  Eq.  (7). 
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To  avoid  build-up  of  turbulence  in  stagnation  regions,  the  model  makes  use  of  a  production  limiter,  as  it  is  possible 
to  find  in  Eq.  (2)  and  Eq.  (3). 
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During  the  upstroke  of  a  pitching  motion,  the  effect  of  an  adverse  pressure  gradient  is  limited,  and  then  it  leads  to  a 
dynamic  stall  angle  that  exceeds  the  experimental  static  stall,  which  increases  the  aerodynamic  coefficients.  However, 
during  the  downstroke,  the  adverse  pressure  gradient  largely  increases,  leading  to  a  reattachment  angle  that  is  smaller 
than  that  of  the  static  case.  As  shown  above,  the  k-co  SST  was  designed  to  capture  recirculation  and  to  accurately  predict 
the  adverse  pressure  gradient  related  to  unsteady  effects  of  transient  flows.  Nevertheless,  since  the  analysis  is  two- 
dimensional,  it  is  possible  that  the  prediction  will  not  be  fully  correct. 


3.  PHYSICAL  MODELLING 


3.1.  Problem  Setting 

The  mesh  used  for  the  analysis  is  a  68500  quad  cell  C-Grid  topology  two-dimensional  mesh.  To  ensure  a  domain- 
independent  solution,  the  size  was  set  to  have  an  upstream  radial  length  of  20  airfoil  chords  from  trailing  edge  and  a 
downstream  length  of  20  chords.  The  domain  is  composed  by  two  sub -domains:  the  external  domain,  with  100  nodes  in 
the  I  direction  and  100  nodes  in  the  J  direction,  with  a  75  nodes  radial  distribution;  the  internal  circular  domain,  with  a  5 
chords  radial  dimension,  has  400  nodes  distributed  in  the  circumferential  direction  and  125  nodes  in  the  radial  direction. 
The  use  of  an  internal  circular  domain  was  to  guarantee  the  mesh  movement  with  respect  to  the  pitching  airfoil  without 
changing  the  quality  of  the  cells  during  motion  in  the  ongoing  simulation. 

The  concept  of  Sliding  Meshes  is  applied  to  the  interfaces  between  the  circular  and  external  domain,  creating  non¬ 
matching  nodes  due  to  the  rotation,  also  known  as  ‘hanging  nodes’.  This  maintains  the  accuracy  of  the  flow  prediction 
close  to  the  wall.  Since  the  nodes  are  supposed  to  have  a  steady  position  in  reference  to  their  moving  frame,  no 
smoothing  dynamic  mesh  method  is  necessary  and  the  quality  does  not  change. 


(a)  (b)  (c) 

Figure  1.  (a)  C-Grid  computational  domain  mesh;  (b)  Internal  circular  domain  mesh;  (c)  Mesh  interface  on 

domain  connection  -  ‘hanging  nodes’. 

In  order  to  avoid  conservation  problems  the  connecting  walls  between  the  domains  are  set  as  interfaces,  so  the  fluid 
will  flow  without  changes  through  it.  This  condition  is  set  to  keep  the  nodes  and  cell  in  the  inner  boundary  of  the 
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external  domain  static,  while  the  nodes  and  cells  in  the  perimeter  of  the  circular  domain  can  slide  following  the  pitching 
airfoil  movement. 

The  governing  equation  of  a  pitch  airfoil  can  be  described  as  shown  in  the  following  equation. 

«  =  CCmln  +  —  2~amin)  (1  -  COS (<Wt))  (9) 


The  unsteady  motion  is  characterized  by  the  reduced  frequency  of  oscillation;  from  where  it  is  possible  to  calculate 
the  sinusoidal  angular  rate,  as  can  be  found  in  Eq.  (9).  Following  the  work  of  McAlister  et  al.  (1978),  the  most 
significant  parameter  in  the  oscillatory  motion  of  an  airfoil  is  its  pitching  frequency  about  its  quarter-chord  axis.  In  his 
work,  it  is  shown  that  the  hysteresis  loop  enlargement  and  the  dynamic  stall  recovery  is  delayed  as  the  reduced 
frequency  is  increased,  deviating  from  the  static  airfoil  values  for  aerodynamic  coefficients  given  as  function  of  the 
incidence  angle. 
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ooc 

2Uoo 
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The  simulation  was  based  on  the  experiments  from  Berton  et  al.  (2002).  For  the  prescribed  motion,  the  minimum 
and  maximum  incidence  angles  were  6  and  18  degrees,  respectively,  the  reduced  frequency  for  the  experiment  was 
given  as  k  =  0.188  at  a  Reynolds  number  of  105,  for  a  1-meter-chord  NACA  0012  airfoil.  To  achieve  the  best  possible 
results,  static  mesh  test  was  made  for  an  incidences  of  0  and  6  degrees,  leading  to  10%  and  5%  deviation  for  both  lift 
and  drag  aerodynamic  coefficients,  respectively,  in  comparison  with  experimental  results  from  Sheldal  et  al.  (1981).  For 
these  settings,  the  mesh  refinement  level  ensures  y-plus  values  that  are  lower  than  one  in  near-wall  region  across  the 
whole  domain.  Results  are  calculated  from  the  first  20  oscillatory  cycles,  whilst  the  analysis  is  performed  from  the 
average  of  the  last  five  simulated  cycles. 


3.2.  Two-dimensional  Simulation 


The  computational  simulation  results  using  the  k-co  SST  model  are  compared  with  the  experiments  from  Berton  et 
al.  (2002)  and  the  k-s  Chien  and  Spalart-Allmaras  model  two-dimensional  simulations  from  Martinat  et  al.  (2008),  for  a 
pitching  axis  located  at  the  center  of  pressure  (one  quarter  chord  aft  the  leading  edge) .  The  aerodynamic  coefficients  of 
lift  and  drag  are  analyzed  as  functions  of  the  incidence  angle  a,  as  shown  below.  The  different  behavior  of  both  lift  and 
drag  coefficients  prediction  can  be  observed  in  Fig.  2.  It  is  found  that  the  k-co  SST  model  presents  a  behavior  that  is 
close  to  the  one  associated  with  the  Spalart-Allmaras  model,  compared  to  the  simulations  from  Martinat  et  al.  (2008). 
The  k-co  SST  model  underestimates  the  lift  coefficient  at  the  upstroke  phase;  however,  it  has  a  less  critical  prediction  for 
lift  and  drag  on  the  downstroke  phase. 

In  the  work  of  Martinat  et  al.  (2008)  it  is  possible  to  evaluate  the  results  from  a  three-dimensional  DDES  k-co  SST 
model.  The  comparison  of  such  results  with  those  from  the  RANS  k-co  SST  model  two-dimensional  simulation  is 
presented  in  Fig.  3.  It  shows  that  the  flow  approximation  is  fairly  close  for  the  upstroke  phase,  indicating  that  the  lift 
and  drag  coefficients  are  not  affected  by  three-dimensional  effects.  This  could  be  also  concluded  from  the  calculated 
results,  since  an  analysis  of  the  last  ten  simulated  cycles  has  shown  that  the  aerodynamics  coefficients  did  not  change 
for  each  incidence  angle.  The  difference  comes  out  during  the  downstroke  phase,  in  which  the  flow  behavior  is  proven 
to  be  three-dimensional  dependent,  as  concluded  in  works  of  McAllister  et  al.  (1978)  and  Martinat  et  al.  (2008). 

For  the  downstroke,  it  can  be  seen  that  the  present  two-dimensional  simulation  presents  a  smaller  deviation  from 
the  experiment,  when  compared  with  the  three-dimensional  simulation  from  Martinat  et  al.  (2008).  However,  both 
numerical  predictions  display  high  oscillation  characteristics  during  the  downstroke  phase,  showing  that  the  model  is 
not  able  to  capture  all  the  circulation  of  an  oscillatory  flow  with  this  scale  of  complexity.  Of  course,  a  flow  of  such 
complexity  of  unsteadiness  is  not  easy  to  model,  as  turbulence  modelling  can  render  misleading  predictions.  Also,  the 
two-dimensional  approximation  can  delay  separation,  giving  results  more  optimistic  than  observed  experimentally. 
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- Experimental  Data  -  Berton  et  al.  (2002) 

- k-Epsilon  Chien  -  Martinat  et  al.  (2008) 

- Spalart  Allmaras  -  Martinat  et  al.  (2008) 

- k-Omega  SST  0.25 


Angle  of  Attack,  a 


- Experimental  Data  -  Berton  etal.  (2002) 

- k-Epsilon  Chien  -  Martinat  et  al.  (2008) 

- Spalart  Allmaras  -  Martinat  et  al.  (2008) 

- k-Omega  SST  0.25 


Figure  2.  Hysteresis  loops  obtained  with  two-dimensional  RANS  k-co  SST  analysis  on  lift  and  drag  coefficients 
compared  with  the  experimental  results  from  Berton  et  al.  (2002)  and  simulations  from  Martinat  et  al.  (2008). 


- Experimental  Data  -  Berton  et  al.  (2002) 

- DDES  k-Omega  SST  -  Martinat  et  al.  (2008) 

- k-Omega  SST  0.25 


Angle  of  Attack,  a 


— Experimental  Data  -  Berton  etal.  (2002) 

- DDES  k-Omega  SST  -  Martinat  et  al.  (2008) 

- k-Omega  SST  0.25 


Figure  3.  Hysteresis  loops  obtained  with  two-dimensional  RANS  k-co  SST  analysis  on  lift  and  drag  coefficients 
compared  with  the  experimental  results  from  Berton  et  al.  (2002)  and  three-dimensional  DDES  k-co  SST 

simulation  from  Martinat  et  al.  (2008). 

The  flow  behavior  of  the  simulation  also  agrees  with  the  studies  presented  by  McCroskey  (1982)  and  details  of 
unsteady  calculation  and  solutions  can  be  found  in  McCroskey  (1973).  During  the  upstroke,  the  flows  remains  attached 
to  the  profile  up  to  12  degrees.  At  a  14.4  degrees  incidence  angle  two  small  recirculation  regions  are  formed,  one  close 
to  the  leading  and  one  close  to  the  trailing  edge;  from  this  point  the  separation  bubble  and  turbulent  separation  point 
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start  moving  towards  each  other  due  to  the  increase  of  the  angle  of  attack.  As  expected,  the  dynamic  stall  for  unsteady 
pitching  airfoils  is  delayed  by  the  influence  of  its  time  dependent  vortices,  whereas  static  stall  would  be  near  to  occur  at 
14  degrees.  At  an  incidence  close  to  17  degrees,  the  trailing  edge  turbulent  separation  and  the  separation  bubble 
generates  a  large  area  of  recirculation,  leading  to  the  detachment  of  the  boundary  layer  at  18  degrees  and  consequent 
drastic  loss  of  lift.  The  flow  is  kept  stalled  during  most  part  of  downstroke,  reattaching  close  to  7  degrees,  as  can  be 
observed  in  the  Lift  Coefficient  vs.  Incidence,  in  Fig.  3,  and  also  in  Fig.  4,  where  the  flow  is  reattached  already  at  7.2 
degrees  during  downstroke. 

In  the  downstroke  phase,  where  the  flow  detaches  from  the  airfoil  surface,  it  is  also  necessary  to  consider  that  the 
lower  surface  of  the  airfoil  is  also  acting,  in  the  other  direction,  producing  down  force,  hence  reducing  the  lift  force  in 
the  phase.  This  occurs  due  to  the  shape  of  the  airfoil.  Since  it  is  symmetric,  both  surfaces  act  as  predicted  by  the  Coanda 
Effect.  The  Coanda  Effect  explains  the  bending  of  fluids  around  an  object  due  to  its  viscosity,  which  makes  a  fluid  thick 
and  makes  it  stick  to  a  surface  (Anderson  &  Eberhardt,  2001).  Due  to  the  difference  of  speed  between  fluid  parcels  of 
the  boundary  layer  and  in  it  vicinity,  shear  forces  are  created,  which  help  attach  the  flow,  and  also  force  it  to  bend  in  the 
direction  of  the  slower  layer,  the  one  close  to  the  wall,  trying  to  wrap  around  the  object.  Hence,  the  Coanda  Effect 
explains  why  there’s  no  lift  on  a  symmetrical  airfoil  at  zero  angle  of  attack,  but  when  it  reaches  a  positive  incidence 
attitude,  the  flow  attaches  to  the  upper  surface  and  due  to  the  shape  of  the  airfoil,  the  flow  is  bent  in  the  trailing  edge, 
generating  down  wash,  which  causes  the  lift  force  over  the  airfoil. 

In  the  downstroke  phase,  the  flow  is  better  attached  to  the  lower  surface  of  the  airfoil,  as  shown  in  Fig.  4,  which 
causes  the  bending  of  the  flow  upwards,  also  generating  a  ‘upwash’  at  the  trailing  edge.  Due  to  vortex  formation  at  the 
trailing  edge  and  stall  condition,  the  boundary  layer  is  detached  from  the  upper  surface,  causing  the  abrupt  loss  of  lift. 
When  the  flow  is  reattached  near  to  7  degrees  downstroke,  the  lower  surface  return  to  bend  the  air  down  in  the  trailing 
edge,  as  expected  for  a  positive  attitude  incidence. 

Even  though  the  k-co  SST  two-dimensional  simulation  overestimated  the  aerodynamic  coefficients  prediction  by 
some  extent,  by  comparing  Fig.  3  it  is  possible  to  see  that  the  prediction  of  the  reattachment  incidence  was  the  closest 
one  to  the  experimental  reference.  Considering  the  two-dimensional  simulations  of  Martinat  et  al.  (2002),  the  latest  k-co 
SST  simulation  matched  accurately  with  the  experiment  downstroke  reattachment  incidence.  For  the  analyzed  model 
performance,  it  is  necessary  to  account  for  the  influence  of  the  SST  limiter,  as  seen  in  Eq.  (8),  which  can  create  a 
reduction  of  the  eddy  viscosity  during  simulation  (Martinat  et  al. ,  2002). 


Figure  4.  Vorticity  for  a  pitching  airfoil  using  k-co  SST  2D  model;  min  =  0.0  [1/s];  max  =  5.0  [1/s]. 
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There  are  many  factors  that  can  affect  the  flow  behavior  under  unsteady  condition.  As  found  in  experiments  by 
McAlister  (1978),  not  only  the  reduced  frequency  have  a  major  importance  on  determining  the  flow  characteristics  and 
vortex  dimensions,  but  also,  other  factors  as  the  thin  airfoil  shape,  flow  velocity  and  rotation  axis.  To  study  the 
influence  of  the  rotation  axis  on  the  vortex  formation  a  test  was  made  in  order  to  compare  results  from  simulations  in 
which  such  axis  was  located  at  the  center  of  pressure  (one  quarter  chord  aft  the  leading  edge),  and  at  the  middle  chord. 

The  influence  of  the  rotation  axis  analysis  was  made  in  order  to  compare  the  results  from  the  k-co  SST  model.  It  is 
shown  by  Fig.  5  that  the  hysteresis  loops  for  lift  and  drag  coefficients  presents  a  close  behavior  for  both  tested  cases. 
Concerning  lift  coefficients  during  pitching  cycle,  when  the  rotation  axis  is  at  the  middle  chord,  a  thinner  hysteresis 
loop  is  observed  in  comparison  with  the  simulations  performed  with  the  rotation  axis  located  at  the  center  of  pressure  of 
the  airfoil.  The  hysteresis  loop  for  the  drag  coefficients  is  larger  when  the  rotation  axis  lies  at  the  middle  chord,  for  high 
incidence  angles. 


- k-Omega  SST  0.25 

- k-Omega  SST  0.50 


- k-Omega  SST  0.25 

- k-Omega  SST  0.50 


Figure  5  -  Hysteresis  loops  obtained  with  two-dimensional  RANS  k-o  SST  computational  analysis  on  lift  and 
drag  coefficients  for  the  rotation  axis  at  the  Center  of  Pressure  (x  =  0.25*chord)  and  at  the  Middle  Chord  (x  = 

0.50*chord). 

Both  lift  and  drag  coefficients  increase  with  an  increase  of  the  velocity  in  the  near- wall  region.  During  the  upstroke 
phase  the  flow  acceleration  is  more  accentuated  for  the  pitching  axis  at  the  center  of  pressure,  which  is  responsible  for 
the  larger  lift  coefficient  for  this  phase.  The  scenario  where  the  rotation  axis  is  at  the  middle  chord  presents  a  lower  lift 
coefficient  curve  in  the  upstroke;  however,  it  is  possible  to  see  a  raise  in  the  drag  coefficient  during  this  phase,  which 
happens  due  to  vortex  formation  in  the  trailing  edge  region.  This  vortex  is  responsible  for  the  formation  of  the  trailing 
edge  turbulent  separation  region,  causing  the  detachment  of  the  boundary  layer  in  that  area.  This  implies  a  decrease  of 
the  boundary  layer  attached  area  on  the  upper  surface,  which  is  responsible  for  the  decrease  in  lift  and  the  increase  in 
drag  generation. 

At  higher  incidence  angles,  the  separation  bubble  that  forms  close  to  the  leading  edge  reveals  itself  more 
accentuated  for  the  center  of  pressure  rotation  axis  condition.  The  boundary  layer  is  attached  to  the  upper  surface  of  the 
airfoil  until  the  incidence  is  greater  than  17  degrees;  at  this  point,  the  bubble  reaches  the  trailing  edge,  causing  a  full 
detachment  of  the  boundary  layer,  hence  leading  to  an  abrupt  lift  loss. 

For  the  middle  chord  case,  the  formation  of  the  separation  bubble  close  to  the  leading  edge  is  delayed.  Unlike  the 
center  of  pressure  rotation  axis  case,  in  this  situation,  in  the  end  of  the  upstroke  phase,  the  bubble  did  not  reach  the 
trailing  edge;  this  will  only  occur  in  the  downstroke  phase  after  reaching  a  17  degrees  angle  of  attack,  causing  a  full 
detachment  of  the  boundary  layer.  Hence,  an  abrupt  lift  loss  occurs  at  an  incidence  angle  close  to  17  degrees  during  the 
downstroke  phase.  The  flow  will  reattach  at  a  low  incidence  angle,  of  about  7  degrees,  as  also  predicted  by  the  center  of 
pressure  pitching  case. 

Both  cases  were  able  to  predict  the  same  reattachment  incidence  angle,  close  to  7  degrees,  as  also  occur  in  the 
experiments  from  Berton  et  al.  (2002).  The  lack  of  experimental  data  for  the  pitching  axis  at  the  middle  chord  prevents 
a  full  analysis  of  the  flow  prediction  for  the  case.  However,  the  simulation  for  the  pitching  axis  at  the  center  of  pressure 
produced  results  that  can  be  accurately  used  for  the  analysis  of  flow  behavior.  Therefore,  it  is  possible  to  assess  first 
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impressions  on  the  influence  of  the  pitching  axis  position  by  comparing  results  from  the  performed  k-co  SST  two- 
dimensional  simulations. 


a  =  17°Downstroke  0.25 


a  =  15 0  Downstroke  0.25 


a  =  12 °Downstroke  0.25 


a  =  7.2° Downstroke  0.25 


a  =  17°Downstroke  0.50  a  =  15° Downstroke  0.50 


a  =  12 °Downstroke  0.50 


a  =  7.2 0 Downstroke  0.50 


Figure  6  -  Velocity  stream  lines  for  a  pitching  airfoil  rotation  at  the  Center  of  Pressure  (x  =  0.25*chord)  and  at 

the  Middle  Chord  (x  =  0.50*chord). 


4.  CONCLUSIONS 

Simulations  were  carried  out  on  a  NACA0012  airfoil  using  the  commercial  software  Fluent  6.4  and  using  the 
RANS  model  k-co  SST  in  order  to  study  the  behavior  of  flow  under  an  unsteady  airfoil  flight  condition.  The  results  were 
compared  with  experimental  data  and  RANS  and  DDES  turbulence  models  to  better  understand  the  behavior  and  check 
accuracy  on  the  flow  characteristics  prediction  of  a  k-co  SST  two-dimensional  simulation  with  Reynolds  Number 
of  10s. 

The  use  of  a  Sliding  Mesh  concept  is  proven  to  be  reliable  for  steady  flows,  as  the  comparison  with  other  results 
cases  from  different  turbulence  models  have  shown  similarity  on  flow  behavior  and  results.  The  mesh  refining  close  to 
the  near-wall  region  is  crucial  for  achieving  an  accurate  modelling,  and  the  use  of  a  sliding  mesh  for  unsteady  flows 
simulation  can  assure  results  as  reliable  as  the  usual  hybrid  tri/quad  dynamic  meshes,  since  it  is  possible  to  maintain  the 
mesh  quality  at  near- wall  regions  for  quad  cells  during  the  mesh  motion. 

Although  the  described  model  was  capable  of  predicting  accurately  well  dynamic  stall  and  the  reattachment 
downstroke  incidence  angle,  the  aerodynamic  coefficients  calculated  for  the  flow  under  a  pitching  condition  deviated 
from  the  experimental  data,  which  occurs  due  to  the  complexity  of  modelling  a  flow  with  this  level  of  unsteadiness  with 
a  two-dimensional,  RANS  approach.  However,  the  behavior  for  the  vorticity  cycles  matched  with  those  expected  from 
literature  results.  Hence,  it  is  possible  to  conclude  that  the  two-dimensional  study  of  a  pitching  airfoil  can  provide  an 
appropriate,  simplified  understanding  of  the  behavior  of  an  unsteady  flow  and  its  vortex  formation  pattern.  The  two- 
dimensional  simulation  was  also  important  to  enlighten  the  Coanda  Effect  during  the  downstroke  phase,  acting  in  the 
lift  loss  at  dynamic  stall. 

The  study  of  the  influence  of  the  pitching  axis  position  in  unsteady  flows  consisted  in  an  initial  step  for  the 
understanding  of  the  variables  that  can  alter  flow  behavior.  For  a  more  deep  understanding,  other  positions  for  the 
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pitching  axis  should  be  analyzed,  rather  than  just  two  cases.  Nevertheless,  the  conducted  preliminary  study  shows  a 
decrease  of  the  dynamic  stall  incidence  angle  and  almost  the  same  reattachment  incidence  for  the  downstroke  phase 
when  one  compares  pitching  movements  about  the  half  chord  and  the  center  of  pressure  locations,  respectively. 

One  also  concludes  that,  even  though  RANS  models  have  been  proved  to  predict  accurately  well  many  industrial 
and  aerodynamic  flows,  it  is  still  hard  to  complex,  unsteady  flow.  The  models  are  able  to  predict  some  of  its 
characteristics,  yet  not  fully.  In  any  event,  the  results  have  been  critically  important  for  the  flow  behavior 
understanding,  where  even  with  the  overprediction  of  some  characteristics,  the  qualitative  behavior  matched  the 
expected  one. 
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Resumo:  Este  texto  apresenta  alguns  dos  principais  estudos  relacionados  a  tecnologia  de  aeronaves  de  geometria  adap- 
tativa.  Inicialmente,  uma  breve  contextualizagao  histdrica  e  apresentada.  Na  sequencia,  relagdes  qualitativas  tipo  causa- 
efeito  sao  consideradas  para  a  constatagao  de  quais  beneficios  podem  resultar  da  adogao  de  geometria  variavel  em 
aeronaves.  For  fim,  sao  considerados  alguns  dos  principais  trabalhos  de  pesquisa  realizados  acerca  do  assunto ,  algumas 
iniciativas  tomadas  pelo  setor  privado,  bem  como  conceitos  relatados  na  literatura.  As  principais  conclusoes  indicam  que 
grande  s  avangos  ja  foram  conquistados,  como  comprovam  implementagoes  funcionais  observadas  em  avioes  de  comb  ate, 
mas  que  esforgos  de  pesquisa  e  desenvolvimento  sao  necessarios  para  que  projetos  de  aeronaves  de  geometria  adaptativa 
tornem-se  eficazes,  seguros  e  pratica  comum,  dada  a  natureza  multidisciplinar  do  tema. 

Palavras-chave:  Aeronaves  de  geometria  adaptativa,  materiais  inteligentes,  estruturas  inteligentes,  estruturas  adaptati- 
vas. 

1.  INTRODUCAO 

Desde  os  primordios  da  humanidade,  o  homem  observa  e  aprende  com  a  natureza.  O  desenvolvimento  do  voo  au- 
tossustentado,  por  exemplo,  comegou  no  final  do  Seculo  XIX,  e  imcio  do  Seculo  XX,  tendo  como  principal  inspiragao 
o  voo  de  passaros.  E  claro  que  a  revolugao  industrial,  bem  como  a  compreensao  de  fenomenos  fisicos  relacionados  a 
mecanica  dos  fluidos,  ajudaram  a  alavancar  o  campo  da  aeronautica  em  momentos  seguintes.  Por  outro  lado,  os  primeiros 
projetos  aeronauticos  foram  fortemente  influenciados  pela  imaginagao  humana.  Prova  disso  sao  desenhos  do  inventor  e 
engenheiro  frances  Clement  Ader  (1841-1925),  que  datam  de  antes  de  1900;  prototipos  chegaram  a  ser  construidos,  mas 
nunca  voaram  [Fig.  (1)].  Um  contemporaneo  de  Ader  foi  o  engenheiro  alemao  Otto  Lilienthal  (1848-1896),  idealizador  de 
planadores;  com  estes,  registrou  mais  de  2000  voos  entre  os  anos  de  1891  a  1896,  quando  de  um  acidente  fatal  [Fig.  (2)]. 

Fica  clara,  se  analisados  estes  e  outros  casos,  a  tentativa  de  algar  voo  por  meio  de  inspiragao  na  natureza.  Animais  que 
detem  a  capacidade  de  voar  o  conseguem  devido  a  evolugao  que  apresentaram  ao  longo  de  milhoes  de  anos.  O  homem, 
por  sua  vez,  tentou  copiar  algumas  das  caracteristicas  mais  marcantes  destes  animais,  como  as  alteragoes  geometricas 
associadas  as  asas  de  aves,  por  exemplo.  Essas  sao  utilizadas  para  levantar  voo,  gerar  propulsao,  garantir  estabilidade, 
entre  outros.  Alem  disso,  alteragoes  geometricas  tambem  permitem  a  uma  ave  alterar  entre  voo  nivelado,  a  realizagao  de 
mergulho  para  caga,  execugao  de  manobras  furtivas,  etc..  Para  exemplificar,  o  falcao  peregrino  e  mostrado  na  Fig.  (3)  em 
situagoes  distintas  de  voo;  esta  ave  e  capaz  de  realizar  sobrevoo  por  longos  periodos  de  tempo,  e  entao  mergulhar  em  alta 
velocidade  para  a  captura  de  alimento. 

Tomando  por  base  o  voo  de  passaros,  os  irmaos  Wright,  ainda  em  1902,  propuseram  o  uso  de  torgao  de  asa  para 
realizagao  do  controle  de  voos  autossustentados.  Uma  patente  foi  concedida  para  este  tipo  de  aplicagao,  e  o  principio  foi 
usado  no  Wright  Flyer,  com  a  adogao  de  um  sistema  de  cabos  e  polias,  acionado  pela  movimentagao  lateral  do  piloto,  que 
se  posicionava  deitado  no  prototipo  de  aviao. 

Com  o  passar  do  tempo,  varias  outras  patentes  que  tinham  como  inspiragao  a  natureza  para  auxiliar  no  controle  de 
aeronaves  foram  registradas.  A  busca  por  maiores  velocidades  de  voo,  contudo,  aliada  a  flexibilidade  das  estruturas 
entao  empregadas  para  construgao  aeronautica,  ocasionou  acidentes  aereos  devidos  a  fenomenos  aeroelasticos,  como 
divergencia  e  flutter.  Assim,  logo  na  sequencia  da  la  Guerra  Mundial,  aeronaves  foram  desenvolvidas  com  materiais  mais 
rigidos,  possibilitando  maiores  velocidades,  mas  limitando  a  aplicabilidade  de  sistemas  inspirados  na  natureza  para  o 


Figura  1 .  Clement  Ader  e 
um  de  seus  prototipos. 


Figura  2.  Otto  Lilienthal  e  Figura  3.  Falcao  peregrino  em  diferentes  configuragoes 
um  de  seus  planadores.  de  voo,  como  se  ve  atraves  de  suas  asas  e  cauda. 
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controle  de  voo,  por  exemplo.  Com  isso,  superficies  auxiliares,  como  flaps,  slats,  freios  aerodinamicos,  etc.,  tornaram-se 
cada  vez  mais  populares  e  presentes  na  aviagao,  desempenhando  papel  de  controle  ou  de  hiper- sustentagao,  conforme  o 
caso  (Weisshaar,  2013). 

Independente  da  forma  como  acontece,  e  evidente  a  necessidade  de  reconfiguragao  geometrica  para  o  voo  de  animais, 
ou  de  aeronaves,  seja  para  possibilitar  controle,  propulsao,  modificagao  de  proposito,  etc..  A  adogao  de  superficies 
auxiliares  em  aeronaves  permanece  ate  hoje  devido  a  seu  aperfeigoamento,  e  a  nao  necessidade  de  sua  substituigao.  Por 
outro  lado,  hoje  ha  a  preocupagao  de  se  produzir  aeronaves  que  sejam  cada  vez  mais  economicas,  ao  menos  na  aviagao 
civil.  Flaps,  por  exemplo,  levam  a  descontinuidades  na  superficie  de  sustentagao,  o  que  implica  maior  esforgo  devido 
ao  arrasto  aerodinamico,  e  a  um  subsequente  maior  gasto  energetico.  Alem  disso,  superficies  auxiliares  tradicionais 
tambem  assumem  responsabilidade  por  grande  parte  do  ruido  gerado  nas  etapas  de  decolagem  e  de  aterrissagem  de  jatos 
comerciais.  Por  estes  motivos,  alternativas  tern  sido  buscadas  para  o  desenvolvimento  de  asas  com  geometria  adaptativa, 
isto  e,  capazes  de  ter  sua  forma  alterada,  como  realizam  aves,  por  exemplo. 

Ao  longo  do  Seculo  XX,  este  tipo  de  abordagem  nao  foi  completamente  esquecida,  como  comprovam  o  contmuo 
registro  de  patentes,  o  desenvolvimento  de  prototipos,  e  a  construgao  de  aeronaves  com  asas  de  enflechamento  variavel, 
por  exemplo  (Weisshaar,  2013).  Contudo,  apesar  de  viabilidade  ter  sido  comprovada  em  varios  casos,  em  outros  o 
acrescimo  de  peso  e/ou  custos  resultavam  proibitivos,  muito  em  fungao  dos  sistemas  empregados  para  atuagao.  Por 
este  motivo,  a  consideragao  de  asas  com  geometria  adaptativa  mostra-se  vantajosa  quase  que  exclusivamente  quando 
requisitos  conflitantes  devem  ser  alcangados  pelo  projeto  aeronautico.  Por  outro  lado,  aeronaves  que  dai  resultam  muito 
provavelmente  serao  superiores  em  relagao  a  outras  alternativas,  geralmente  por  serem  versateis  gragas  a  possibilidade  de 
reconfiguragao. 

Outro  ponto  que  merece  destaque  para  o  aumento  de  interesse  relacionado  a  aeronaves  com  geometria  adaptativa  e 
o  desenvolvimento  dos  denominados  materiais  inteligentes  (Leo,  2007).  Estes  podem  ser  empregados  como  atuadores 
que  dispoem  de  vantagens  por  serem  compactos,  possibilitar  atuagao  distribuida,  desempenhar  simultaneamente  fungao 
estrutural,  e  possuirem  alta  densidade  energetica  ou  altas  taxas  de  resposta,  conforme  o  caso.  Sistemas  de  atuagao  que 
empregam  tais  materiais  podem  levar  a  uma  redugao  de  peso  se  considerados  para  comparagao  aqueles  sistemas  tradici¬ 
onais  de  natureza  hidraulica  e/ou  pneumatica.  Exemplos  de  materiais  inteligentes  compreendem  materiais  piezeletricos, 
ligas  com  memoria  de  forma,  fluidos  magneto-reologicos  e  eletro-reologicos,  etc.. 

Considerando  o  exposto  anteriormente,  o  presente  resumo  expandido  busca  apresentar  alguns  dos  principais  avangos 
relacionados  a  tecnologia  envolvida  com  aeronaves  de  geometria  adaptativa.  Na  sequencia,  enfoque  e  dado  as  possibi- 
lidades  ate  entao  consideradas  para  modificagao  geometrica.  Projetos  de  pesquisa  e  conceitos  expostos  na  literatura  sao 
entao  brevemente  considerados.  Por  fim,  conclusoes  sao  apresentadas. 

2.  O  QUE  ADAPTAR  E  POR  QUE? 

Como  apontado  por  Jha  e  Kudva  (2004),  varias  sao  as  opgoes  para  reconfiguragao  geometrica  de  asas  de  aeronaves. 
De  acordo  com  as  teorias  classicas  empregadas  em  seu  projeto  [como  detalha  Raymer  (2012),  por  exemplo],  influencias 
podem  ser  estabelecidas  de  modo  qualitativo  como  apresentado  na  Tab.  (1).  Por  outro  lado,  o  projeto  de  aeronaves  com 
geometria  adaptativa  deve  ser  conduzido  com  cautela,  tendo  em  vista  que,  como  ja  mencionado,  pode  haver  acrescimo  de 
peso,  e  reconfiguragao  pode  sobrecarregar  um  ou  outro  sub-sistema  com  demandas  adicionais. 


Tabela  1.  Efeitos  oriundos  de  alteragoes  geometricas  na  asa  sobre  o  desempenho  de  aeronaves  (Jha  e  Kudva,  2004). 


Parametro 

Variagao 

Efeitos  resultantes  (com  parametros  restantes  mantidos  fixos) 

Area  planar  da  asa 

+ 

Aumento  na  sustentagao  e  na  capacidade  de  carga 

Diminuigao  do  arrasto  parasita 

Razao  de  aspecto 

+ 

Aumento  de  L/D,  do  tempo  de  loiter,  na  distancia  de  cruzeiro  e  nas  taxas  de  curvas;  Diminuigao 
nos  requisitos  para  motores 

Aumento  da  velocidade  maxima;  Diminuigao  no  arrasto  parasita 

Angulo  de  diedro 

+ 

Aumento  do  momento  de  rolagem  e  da  estabilidade  lateral 

Aumento  da  velocidade  maxima 

Enflechamento 

+ 

Aumento  do  Macritico  e  do  efeito  diedro;  Diminuigao  do  arrasto  devido  a  ondas  de  choque 
Aumento  do  (Cl) max 

Razao  de  afilamento 

± 

Eficiencia  da  asa  (distribuigao  da  sustentagao)  e  arrasto  induzido 

Distribuigao  de  torgao 

± 

Previne  o  stall  de  ponta  de  asa  e  altera  a  distribuigao  de  sustentagao 

Curvatura  do  aerofolio 

+ 

Modifica  o  angulo  de  sustentagao  nula,  a  eficiencia  do  aerofolio  e  o  comportamento  de  separagao 

Razao  espessura/corda 

+ 

Aerofolio  com  melhor  desempenho  em  baixas  velocidades 

Aerofolio  com  melhor  desempenho  em  altas  velocidades 

Raio  do  bordo  de  ataque 

+ 

Aerofolio  com  melhor  desempenho  em  baixas  velocidades 

Aerofolio  com  melhor  desempenho  em  altas  velocidades 

Distribuigao  de  espessura 

+ 

Modifica  caracterfsticas  do  aerofolio,  e  a  transigao  laminar/turbulento 
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(a)  Mission  Adaptive  Wing 


(b)  Smart  Wing  Program 


(c)  Aircraft  Morphing  Program 


(d)  Morphing  Aircraft  Structures  Program  (e)  Asa  oblfqua 


Figura  4.  Principals  iniciativas  governamentais  dos  EUA  relacionadas  a  aeronaves  de  geometria  adaptativa. 


Das  possibilidades  apontadas  na  Tab.  (1),  umas  sao,  por  razoes  construtivas,  de  implementagao  claramente  mais 
facil  que  outras.  Barbarino  et  al.  (2011)  confirmam,  por  exemplo,  que  a  maioria  dos  estudos  recentes  que  lidam  com 
aeronaves  de  geometria  adaptativa  visam  alteragoes  na  curvatura  de  aerofolios  e  na  distribuigao  de  torgao  ao  longo  da 
asa.  Modificagoes  de  enflechamento  e  de  envergadura  (alteragoes  simultaneas  da  area  planar  da  asa  e  de  sua  razao  de 
aspecto)  ja  foram  relativamente  bastante  estudadas,  e  se  encontram  num  estado  maduro,  a  ponto  de  ja  terem  permitido 
implementagoes  funcionais.  Poucos  estudos  interessaram-se  por  variagoes  no  angulo  de  diedro  e  na  flexao  ao  longo  da 
envergadura  de  asas.  Alteragoes  na  corda  tern  sido  menos  investigadas,  ja  que  trazem  poucos  beneficios  e  sao  muito 
desafiadoras  do  ponto  de  vista  estrutural. 


3.  TRABALHOS  RECENTES 

Nos  ultimos  anos,  incontaveis  trabalhos  cientificos  foram  publicados  acerca  de  aeronaves  com  geometria  adaptativa. 
Prova  disto  e  a  publicagao  de  varias  revisoes  bibliograficas,  assim  como  de  um  livro  a  respeito  do  tema  (Weisshaar,  2006; 
Thill  et  al .,  2008;  McGowan  et  al .,  2009;  Sofia  et  al .,  2010;  Gomez  e  Garcia,  2011;  Barbarino  et  al ,  2011;  Taha  et  al ., 
2012;  Valasek,  2012;  Vasista  et  al .,  2012;  Weisshaar,  2013;  Sinapius  et  al .,  2014;  Barbarino  et  al .,  2014). 

Grandes  projetos  de  pesquisa  realizados  nos  Estados  Unidos  da  America  incluem:  o  programa  Mission  Adaptive  Wing 
(1975-1988);  o  Smart  Wing  Program ;  o  Aircraft  Morphing  Project  (NASA  Langley,  1994-2004);  e  o  Morphing  Aircraft 
Structures  Program  (DARPA,  2002-2007).  La,  ha  ainda  a  investigagao  do  uso  de  asas  obliquas,  que  conta  com  parceria 
entre  a  NASA  e  a  DARPA,  entre  outras  iniciativas.  A  Fig.  (4)  ilustra  parte  destas  iniciativas. 

Na  Europa,  o  Centro  Aeroespacial  Alemao  (DLR)  e  o  principal  responsavel  por  diversos  conceitos,  como  o  finger  e 
o  belt-rib ,  mostrados  na  Fig.  (5).  Tambem  foram  pesquisados  pelo  instituto  arqueamento  variavel  no  bordo  de  fuga,  e  a 
utilizagao  de  um  bump  para  redugao  de  arrasto  em  regime  trans-sonico.  Outros  estudos  tiveram  como  enfoque  um  bordo 
de  ataque  adaptativo;  um  prototipo  fora  construido  para  validagao  bem-sucedida  do  conceito. 

Empresas  privadas  tambem  mostram  certo  interesse  no  desenvolvimento  de  aeronaves  adaptativas.  A  FlexSys  desen- 
volveu  o  adaptive  compliant  trailing  edge  a  partir  de  mecanismos  flexiveis;  as  ultimas  noticias  a  respeito  do  prototipo 
de  produto  dao  conta  da  realizagao  de  testes  em  voo.  A  Festo  desenvolveu  o  SmartBird ,  que  consiste  em  uma  aeronave 
nao-tripulada  controlada  remotamente,  e  que  teve  seus  mecanismos  de  funcionamento  e  de  voo  completamente  baseados 
na  mecanica  exibida  por  passaros.  Uma  parceria  entre  a  Bell  e  a  Boeing  resultou  em  uma  aeronave  extremamente  versatil, 
o  V-22  Osprey,  que  funciona  como  pega  estrategica  para  as  forgas  armadas  norte-americanas. 

Outros  varios  conceitos  tambem  foram  apresentados  na  literatura.  Estes  compreendem  a  adogao  de  fios  de  liga  com 
memoria  de  forma;  a  aplicagao  de  motores  ultrassonicos  piezeletricos;  o  emprego  de  macro-fiber  composites ;  o  uso  de 
materiais  compostos  bi-  ou  multi-estaveis;  entre  outros.  Varios  trabalhos  tambem  se  dedicam  a  modelagem  e  a  metodos 
para  analise  de  sistemas  aeronauticos  com  geometria  adaptativa,  tendo  em  vista  a  necessidade  de  projeto  adequado. 

4.  CONCLUSAO 

Os  varios  trabalhos  de  pesquisa  divulgados  na  literatura  cientifica  comprovam  a  relevancia  dos  estudos  acerca  de 
aeronaves  de  geometria  adaptativa.  Varios  avangos  ja  foram  conseguidos,  como  comprovam  implementagoes  operacionais 
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Figura  5.  Conceitos  desenvolvidos  pelo  DLR. 


Figura  6.  Envoi vimento  da  iniciativa  privada 
com  aeronaves  de  geometria  adaptativa. 


que  envolvem,  por  exemplo,  enflechamento  variavel  para  operagao  de  cagas  em  regimes  sub-,  trans-  e  super- sonico. 
Pesquisas  mais  recentes  buscam  a  integragao  de  alternativas  para  atuagao,  visando  redugao  de  peso,  e  projetos  mais 
adequados,  pela  distribuigao  de  esforgos  de  atuagao. 

Por  outro  lado,  aeronaves  com  geometria  variavel  ainda  devem  levar  certo  tempo  para  se  tornar  uma  opgao  viavel. 
Barbarino  et  al.  (2011),  por  exemplo,  mencionam  que,  de  acordo  com  a  NASA,  20  a  30  anos  serao  necessarios  para  que 
as  aeronaves  tornem-se  mais  parecidas  com  passaros,  e  deixem  de  ter  superficies  de  controle  discretas  para  alterar  sua 
forma  de  maneira  mais  fluida  e  continua.  Ademais,  a  natureza  multidisciplinar  associada  a  estruturas  inteligentes  adap- 
tativas  acarreta  aumento  significativo  de  complexidade  ao  projeto  aeronautico.  Assim,  varios  avangos  ainda  se  mostram 
necessarios  para  a  adogao  de  asas  com  curvatura  variavel,  ou  com  distribuigao  de  torgao  adaptativa,  por  exemplo. 
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